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Total  Evaporative  Resistance  of  Selected  Clothing  Ensembles 

Victor  Caravello 

ABSTRACT 

With  regard  to  heat  stress,  the  limiting  factor  inherent  in  clothing  ensembles  is  the 

total  evaporative  resistance.  For  the  same  work  demands,  clothing  with  higher 

evaporative  resistance  impedes  the  ability  to  cool  by  sweat  evaporation.  Knowing  the 

evaporative  resistance  provides  a  means  to  compare  candidate  ensembles.  Further,  a 

value  for  evaporative  resistance  means  that  a  rational  method  can  be  used  to  assess  the 

heat  stress  exposure.  Evaporative  resistance  of  five  clothing  ensembles  (cotton  work 

clothes,  cotton  coveralls,  and  three  coveralls  of  particle  barrier,  liquid  barrier  and  vapor 

barrier  properties)  was  determined  empirically  from  wear  tests  during  two  study  phases. 

For  Phase  1,  the  metabolic  rate  was  held  constant  at  160  W/m2,  and  three  levels  of 

humidity  (20,  50,  70%  rh)  were  explored.  Fourteen  heat-acclimated  participants  (9  men 

and  5  women)  completed  trials  for  all  combinations  of  clothing  ensemble  and 

environment.  In  the  Phase  2  study,  the  humidity  was  held  constant  at  50%  rh,  and  three 

levels  of  metabolic  rate  (1 14,  176,  250  W/m2)  were  explored.  Fifteen  heat-acclimated 

participants  (11  men  and  4  women)  completed  trials  for  all  combinations  of  clothing 

ensemble  and  environment.  The  data  from  both  phases  were  analyzed  separately  using 

« 

ANOVA.  Significant  differences  were  found  among  ensembles  (p<0.0001).  The  vapor 
barrier  ensemble  had  the  highest  resistance  at  0.026  kPa  m2/W.  The  liquid  barrier  was 


v 


next  at  0.01 8;  followed  by  the  particle  barrier  and  cotton  coveralls  at  0.016.  Work 
clothes  was  0.014  kPa  m2/W.  Pair-wise  comparisons  adjusted  for  multiple  comparisons 
were  used  to  locate  the  differences  among  ensembles.  Vapor  and  liquid  barrier  ensembles 
were  found  to  be  significantly  different  from  other  ensembles.  Data  from  both  studies 
support  the  conclusion  that  there  are  differences  in  evaporative  resistances  among 
selected  ensembles  tested.  From  the  Phase  2  study,  Ensembles  B  -  E  evaporative 
resistances  decreased  from  0.0024  to  0.0094  kPa  m2/W  with  increased  activity.  The 
decreased  evaporative  resistances  in  Phase  2  can  be  explained  by  the  pumping  action 
associated  with  increased  work.  The  relationship  of  Re,T  to  the  difference  of  Pair  -  Pskin 
(AP)  was  explored  and  found  Rej  does  not  remain  constant.  Environment  appeared  to 
have  greater  influence  on  this  relationship. 


vi 


Total  Evaporative  Resistance  of  Selected  Clothing  Ensembles 


Victor  Caravello 

ABSTRACT 

With  regard  to  heat  stress,  the  limiting  factor  inherent  in  clothing  ensembles  is  the 
total  evaporative  resistance.  For  the  same  work  demands,  clothing  with  higher 
evaporative  resistance  impedes  the  ability  to  cool  by  sweat  evaporation.  Knowing  the 
evaporative  resistance  provides  a  means  to  compare  candidate  ensembles.  Further,  a 
value  for  evaporative  resistance  means  that  a  rational  method  can  be  used  to  assess  the 
heat  stress  exposure.  Evaporative  resistance  of  five  clothing  ensembles  (cotton  work 
clothes,  cotton  coveralls,  and  three  coveralls  of  particle  barrier,  liquid  barrier  and  vapor 
barrier  properties)  was  determined  empirically  from  wear  tests  during  two  study  phases. 
For  Phase  1,  the  metabolic  rate  was  held  constant  at  160  W/m2,  and  three  levels  of 
humidity  (20,  50, 70%  rh)  were  explored.  Fourteen  heat-acclimated  participants  (9  men 
and  5  women)  completed  trials  for  all  combinations  of  clothing  ensemble  and 
environment.  In  the  Phase  2  study,  the  humidity  was  held  constant  at  50%  rh,  and  three 
levels  of  metabolic  rate  (1 14,  176, 250  W/m2)  were  explored.  Fifteen  heat-acclimated 
participants  (1 1  men  and  4  women)  completed  trials  for  all  combinations  of  clothing 
ensemble  and  environment.  The  data  from  both  phases  were  analyzed  separately  using 
ANOVA.  Significant  differences  were  found  among  ensembles  (p<0.0001).  The  vapor 
barrier  ensemble  had  the  highest  resistance  at  0.026  kPa  m2/W.  The  liquid  barrier  was 


v 


next  at  0.018;  followed  by  the  particle  barrier  and  cotton  coveralls  at  0.016.  Work 
clothes  was  0.014  kPa  m2/W.  Pair-wise  comparisons  adjusted  for  multiple  comparisons 
were  used  to  locate  the  differences  among  ensembles.  Vapor  and  liquid  barrier  ensembles 
were  found  to  be  significantly  different  from  other  ensembles.  Data  from  both  studies 
support  the  conclusion  that  there  are  differences  in  evaporative  resistances  among 
selected  ensembles  tested.  From  the  Phase  2  study,  Ensembles  B  -  E  evaporative 
resistances  decreased  from  0.0024  to  0.0094  kPa  m2/W  with  increased  activity.  The 
decreased  evaporative  resistances  in  Phase  2  can  be  explained  by  the  pumping  action 
associated  with  increased  work.  The  relationship  of  Rej  to  the  difference  of  Pair  -  Pskin 
(AP)  was  explored  and  found  Re,T  does  not  remain  constant.  Environment  appeared  to 
have  greater  influence  on  this  relationship. 


vi 


TABLE  OF  CONTENTS 


LIST  OF  TABLES . iii 

LIST  OF  FIGURES . iv 

ABSTRACT . v 

INTRODUCTION . 1 

Heat  Exchange . . . 1 

Role  of  Clothing  in  Heat  Balance . 5 

Insulation . 5 

Permeability . ...6 

Ventilation . . . 6 

LITERATURE  REVIEW . 7 

Components  of  Insulation  and  Evaporative  Resistance . 9 

Laboratory  Test  Methods . 12 

Hot  Plate  Method . 12 

Copper  Manikin . . 1 3 

Human  Tests  for  Clothing  Thermal  Characteristics . . . 1 5 

Thermal  Resistance  Values  for  Various  Work  Ensembles . 23 

Hypothesis . 26 

METHODS . 28 

Participants . . 28 

Clothing  Ensembles . . . 29 

Protocols . 30 

Trials . 31 

Critical  Conditions . 34 

Statistical  Analysis . 36 

RESULTS . 40 

Experimental  Data . 40 

Total  Insulation . 41 

Phase  1 . 42 

Phase  2 . 47 

i 


DISCUSSION 


Internal  Validity . . . . . . 

Comparison  to  Other  Studies . 

Total  Insulation . . 

Total  Evaporative  Resistance. 

Phase  1 . . . . . 

Phase  2 . 


CONCLUSION . 

REFERENCES . 

APPENDIX  A:  PARTICIPANT  DATA . 

APPENDIX  B:  EXPERIMENTAL  DATA  -  PHASE  1 . 

APPENDIX  C:  EXPERIMENTAL  DATA  -  PHASE  2 . . . 

APPENDIX  D;  SAS  CODE  AND  ANALYSIS  -  PHASE  1 .......... 

APPENDIX  E:  SAS  CODE  AND  ANALYSIS  -  PHASE  2 .......... 

APPENDIX  F.  JMP  IN  DATA  ANALYSIS  -  PROTOCOLS. . . 

APPENDIX  G.  JMP  IN  DATA  ANALYSIS  -  M2R5 . . 

APPENDIX  H.  GRAPHS  OF  R^T  VERSUS  AP  BY  ENSEMBLE 
ABOUT  THE  AUTHOR . 


. 53 

. . .53 

. 56 

. ...56 

. 58 

...........59 

. 62 

. 64 

. 66 

. 70 

. 72 

. . 81 

. . ..91 

..........116 

. . 132 

. . 143 


. 149 

End  Page 


ii 


LIST  OF  TABLES 


Table  1 .  Comparison  of  Experimentally  Determined  Thermal  Resistance 

Values . 23 

Table  2.  Summary  of  Participant  Characteristics . 29 

Table  3.  Total  Insulation  Values  for  Ensembles . 41 

Table  4.  Phase  1  -  Mean  Rej  Values  with  Standard  Deviations . 42 

Table  5.  Phase  2  -  Mean  R^j  Values  with  Standard  Deviations . 47 

Table  6.  Phase  2  -  Average  Metabolic  Rates . . 50 

Table  7.  Statistical  Analysis  of  M2R5  Dataset . 55 

Table  8.  It  Values  from  Different  Studies . 57 

Table  9.  Rej  Values  from  Different  Studies . 59 

Table  1 0.  Regression  Analysis  -  AP  by  Rej . 62 


iii 


LIST  OF  FIGURES 


Figure  1. 
Figure  2. 
Figure  3. 
Figure  4. 
Figure  5, 
Figure  6. 
Figure  7. 
Figure  8. 
Figure  9. 
Figure  10 
Figure  11, 
Figure  12 
Figure  13 


Typical  Time  Course  For  Tre  During  an  Inflection  Point  Protocol, 

Phase  1  -  Mean  Rej  by  Ensemble . 

Phase  1  -  Mean  Rej-  by  Ensemble  w/o  Ensemble  E . 

Phase  1 — Mean  R^t  Values  by  Ensemble  and  Environment . . 

Phase  1 — Mean  Re  t  Values  by  Environment . . . 

Phase  2  -  Mean  Re/rby  Ensemble.... . ...................................... 

Phase  2  -  Mean  Rejby  Metabolic  Rate . . . . 

Phase  2  -  Average  Metabolic  Rates  by  Ensemble, . . 

Phase  2  —  Average  Re,Tby  Metabolic  Rate  and  Ensemble . 

Comparison  of  M2R5  Mean  Rej. . . 

Comparison  of  M2R5  Mean  Metabolic  Demands . 

Effect  of  Environment:  Re/r  vs.  AP  -  Ensemble  A . 

Effect  of  Metabolic  Rate  on  Rej . 


.18 

.44 

,44 

,45 

.46 

.49 

.49 

.51 

,51 

.54 

.55 

.61 

.63 


iv 


INTRODUCTION 


Personal  protective  clothing  has  become  commonplace  in  many  industries  today. 
While  protective  clothing  provides  protection  from  exposure  to  chemical  and  physical 
agents,  it  may  lead  to  another  health  issue  -  excessive  heat  strain.  Heat  strain,  the 
physiological  adjustment  to  heat  stress,  is  driven  by  work  demand,  environmental  factors 
(such  as  ambient  temperature,  relative  humidity  and  air  movement),  and  clothing 
requirements.  Heat  strain  is  marked  by  increased  body  temperature,  heart  rate  and 
sweating.  Heat  stress  has  been  studied  extensively,  and  one  of  the  critical  factors  that  ties 
protective  clothing  to  heat  stress  is  evaporative  resistance. 


Heat  Exchange 

To  better  understand  the  role  of  clothing  in  heat  stress,  the  workplace  factors 
discussed  above  can  be  described  through  a  thermal  balance  model.  This  model  balances 
net  heat  gained  by  the  body  with  the  required  heat  loss  to  prevent  excessive  heat  buildup; 
that  is,  to  maintain  thermal  equilibrium.  Thermal  balance  is  frequently  described  by 
some  variation  of  Equation  1  [1,2]. 
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S  =  (M-W)  +  R  +  C-E 


(1) 


The  heat  storage  rate  (S)  represents  the  net  heat  gain  to  or  loss  from  the  body.  By 
convention,  body  temperature  increases  when  S  is  positive  and  decreases  when  heat  is 
lost  (S  is  negative).  When  S  =  0,  the  body  is  considered  to  be  in  thermal  equilibrium. 
Heat  is  generated  internally  by  metabolism  (M),  The  rate  of  metabolic  heat  gain  is 
determined  by  the  rate  and  type  of  external  work  performed  by  the  body.  The  total  heat 
generated  by  metabolic  demands  from  the  work  is  equal  to  metabolic  rate  less  the  rate  of 
external  work  performed  (W).  The  rate  of  radiant  heat  transfer  (R)  between  the  skin  and 
the  environment  and  the  rate  of  heat  transfer  between  the  air  and  skin  surface  (C) 
collectively  characterize  dry  heat  exchange.  Positive  values  for  R  +  C  are  a  heat  gain 
while  negative  values  indicate  a  heat  loss.  The  term  E  represents  the  rate  of  evaporative 
cooling  due  to  the  evaporation  of  sweat,  which  is  the  primary  mechanism  for  cooling  the 
body. 


Heat  production  is  determined  by  the  amount  of  metabolic  activity.  At  rest,  the 

body  generates  heat  from  the  energy  produced  to  maintain  basic  body  functions  such  as 

respiration  and  heart  rate  to  supply  the  needed  oxygen  and  nutrients  to  the  cells. 

Metabolic  activity  rises  as  one  becomes  more  active.  This  rise  results  with  higher 

demand  for  oxygen  and  nutrients  accompanied  by  increased  metabolism  at  the  active 

muscles.  With  increased  metabolism,  there  is  increased  heat  production  in  the  muscle. 

The  greater  the  demand,  the  more  internal  heat  is  generated.  With  the  understanding  that 
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R  +  C  have  a  lesser  effect  of  increasing  internal  heat  production,  the  fundamental  link 
between  metabolic  rate  (work  demands)  to  heat  storage  becomes  clear  [3]. 

The  minimal  effect  of  dry  heat  exchange  (R  +  C)  as  compared  to  the  other  terms 
in  Equation  1  becomes  evident  as  these  terms  are  examined  closer.  Looking  first  at  heat 
transfer  rate  by  convection  (C),  as  ambient  air  temperature  is  raised  above  skin 
temperature,  the  rate  of  heat  gain  by  convection  is  increased.  Simply  stated,  C  is  the 
difference  of  the  ambient  air  temperature  and  the  average  skin  temperature  and  modified 
by  the  rate  of  air  movement  over  the  skin.  If  the  ambient  air  is  cooler  than  the  skin,  then 
heat  flows  away  from  the  body.  The  rate  and  direction  of  convective  heat  exchange 
depends  on  the  temperature  gradient  between  the  air  and  the  skin.  The  rate,  but  not 
direction,  is  also  influenced  by  air  motion  and  clothing.  Generally  the  higher  the  air 
motion  or  velocity,  the  greater  the  rate  of  heat  transfer.  Clothing  provides  a  barrier  to  the 
heat  transfer  through  its  insulation,  so  the  more  skin  that  is  covered  and/or  the  thicker  the 
clothing,  the  lower  the  rate  of  convective  heat  transfer. 

The  temperature  of  surrounding  objects  affects  the  radiant  heat  exchange  between 

the  environment  and  the  body.  Surfaces  of  different  temperatures  have  a  net  heat  flow 

from  the  hotter  to  the  cooler  surface  by  thermal  radiation.  The  rate  of  heat  transfer  by 

radiant  heat  (R)  depends  on  two  factors.  The  first  is  the  temperature  gradient  between  the 

skin  and  surrounding  objects.  If  the  average  temperature  of  surrounding  objects  is 

greater  than  skin  temperature,  there  is  a  heat  gain.  Conversely,  if  the  surroundings  have  a 
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lower  average  temperature  than  the  skin,  a  heat  loss  occurs.  The  rate  of  heat  transfer  is 
proportional  to  the  temperature  gradient.  The  second  factor  is  clothing  insulation.  As 
with  convection,  the  thicker  the  clothing  and/or  the  more  skin  that  is  covered,  the  lower 
the  rate  of  heat  transfer.  Also,  if  the  clothing  has  a  reflective  surface,  the  thermal 
radiation  (or  heat)  is  reflected  away.  Once  the  environmental  temperature  exceeds  35°F 
the  body  can  dissipate  heat  only  by  evaporation  [2]. 

Evaporation  of  sweat  from  the  skin  is  the  primary  mechanism  for  losing  excess 
body  heat  during  activity.  However,  there  is  a  limit  to  the  amount  of  evaporative  cooling 
that  can  occur.  This  limit  is  due  to  two  factors;  a  physiological  and  a  physical  limit.  The 
physiological  limit  is  the  amount  of  sweat  that  can  be  produced  over  time.  The  physical 
limit  is  the  maximum  rate  of  evaporative  cooling  (Emax)  that  can  occur.  Emax  is  limited  by 
three  primary  factors.  First,  evaporation  can  only  occur  if  the  water  vapor  pressure  of  the 
skin  (PSk)  is  higher  than  the  water  vapor  pressure  of  the  ambient  air  (Pa).  Humidity  is  the 
ambient  water  vapor  pressure.  As  humidity  increases,  this  gradient  from  skin  to  air  is 
reduced,  and  the  rate  of  evaporative  heat  loss  is  decreased.  The  second  factor  is  air 
movement.  As  air  velocity  increases,  the  boundary  layer  between  the  person  and 
environment  decreases  allowing  for  an  increase  in  evaporation  of  sweat  (Emax  increases). 
The  third  factor  is  clothing.  All  clothing  acts  as  a  barrier  to  evaporation.  The  physical 
characteristic  known  as  the  water  vapor  permeability  of  the  clothing  is  directly 
proportional  to  the  ability  to  evaporate  sweat  (Emax).  Therefore,  as  the  water  vapor 
permeability  decreases,  so  does  Emax- 
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Air  velocity,  generated  by  body  movements  and  air  movement,  is  important  in 
heat  exchange  between  the  body  and  the  environment  because  of  its  role  in  convective 
and  evaporative  heat  transfer.  Increasing  the  air  velocity  can  increase  the  convective  and 
evaporative  heat  exchange  by  forcing  air  between  the  clothing  and  skin. 

Role  of  Clothing  in  Heat  Balance 

Clothing  impedes  heat  exchange  between  the  body  and  the  environment  by 
limiting  dry  heat  exchange  and  evaporative  cooling.  These  effects  can  be  described  in 
further  detail  by  looking  at  three  characteristics  associated  with  clothing:  insulation, 
permeability  and  ventilation  [4], 

Insulation 

Insulation  describes  the  resistance  to  heat  flow  by  convection  and  radiation.  With 
the  environmental  conditions  being  constant,  the  gradient  between  the  skin  and  air 
remains  the  same,  but  the  as  the  insulation  for  an  ensemble  increases,  the  heat  flow  due  to 
radiation  and  convection  decreases.  In  other  words,  dry  heat  exchange  through  the 
clothing  decreases  with  increasing  insulation. 
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Permeability 

Permeability  is  the  ability  of  water  vapor  to  move  through  clothing.  It  affects  the 
amount  of  evaporative  cooling  that  can  occur.  Clothing  with  low  permeability  indicates 
that  evaporation  of  sweat  through  the  clothing  is  reduced,  resulting  in  a  decrease  in 
evaporative  cooling.  Protective  clothing  ensembles  can  vary  over  a  range  from  easily 
permeable  to  essentially  impermeable. 

Ventilation 

Ventilation  occurs  as  ambient  air  moves  through  the  fabric  and/or  through 
clothing  openings  (cuffs,  fasteners,  and  collar).  Clothing  that  allows  air  movement 
increases  convective  and  evaporative  cooling.  Conversely,  if  the  clothing  is  designed  to 
limit  the  movement  of  air  by  being  encapsulating  or  tight  fitting  with  elastic  cuffs,  the 
convective  and  evaporative  cooling  are  limited. 

Although  protective  clothing  ensembles  are  worn  to  protect  workers  from 
biological,  chemical  or  physical  hazards,  the  barrier  posses  another  hazard  to  workers  by 
reducing  the  wearer’s  ability  to  dissipate  internally  generated  heat  through  sweat 
evaporation.  Depending  on  the  environment  and  work  demands,  an  excessive  level  of 
heat  stress  may  result.  Heat  stress  may  cause  reduced  performance  and  increased  risk  of 
accidents  and  heat  injury.  It  is  imperative  to  understand  the  clothing  characteristics,  with 
the  evaporative  resistance  being  the  most  important,  in  order  to  effective  manage  the 

risks  associated  with  wearing  the  protective  clothing. 

6 


LITERATURE  REVIEW 


Protective  clothing  and  environmental  conditions  influence  the  level  of  heat  stress 
a  worker  may  experience.  Understanding  fabric  properties  may  help  predict  how  the 
environment  will  affect  heat  transfer  for  a  selected  ensemble.  Havenith  points  out  the 
importance  of  heat  balance  when  wearing  protective  clothing  [4].  The  goal  is  to  maintain 
the  body  at  around  37°C  by  transferring  excess  heat  from  the  body  to  the  environment. 
Heat  is  produced  through  metabolic  activity  and  protective  clothing  may  hinder  the  loss 
of  the  heat  gained.  Some  important  factors  that  affect  heat  transfer  from  the  body  to  the 
environment  include  the  temperature  (air,  surface,  radiant),  humidity,  wind,  movement, 
and  clothing  insulation  [4].  While  all  of  these  factors  can  affect  heat  transfer,  the  primary 
mechanism  the  body  uses  to  dissipate  heat  is  evaporative  cooling.  Therefore  it  is 
important  to  understand  the  potential  barrier  an  ensemble  may  pose  to  evaporation  of 
sweat.  As  sweat  is  secreted  onto  the  skin,  it  should  evaporate  and  cool  the  body.  The 
rate  of  evaporation  depends  on  the  difference  between  the  water  vapor  pressure  of  the 
skin  and  the  ambient  air  water  vapor  pressure  as  well  as  the  barrier  provided  by  the 
clothing.  This  barrier  interferes  with  the  ability  of  water  vapor  to  pass  from  the  skin 
through  the  ensemble  and  into  the  ambient  air.  Therefore  it  is  important  to  be  able  to 
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distinguish  between  clothing  ensembles  in  terms  of  their  permeability  to  water  vapor. 
Permeability  is  alternatively  expressed  as  total  evaporative  resistance  (Rej). 

In  addition  to  evaporative  cooling,  heat  loss  from  the  skin  to  the  ambient  air 
occurs  by  radiation  and  convection  (dry  heat  exchange).  Dry  heat  exchange  occurs 
because  of  the  temperature  difference  between  the  skin  and  surrounding  air.  As  with 
evaporation  of  water,  dry  heat  also  must  leave  the  skin  and  be  transported  into  and  out  of 
the  clothing  before  the  heat  loss  is  complete.  Therefore,  clothing  may  interfere  with  dry 
heat  exchange.  This  characteristic  of  clothing  is  referred  to  as  insulation. 

The  total  clothing  insulation  (It)  and  the  total  evaporative  resistance  (Rej)  are 
important  characteristics  to  consider  when  comparing  clothing  ensembles.  It  is  an 
attribute  that  accounts  for  a  decrease  in  heat  flow  due  to  total  insulation  provided  by  the 
clothing  and  the  air  layer  between  the  skin  and  clothing.  The  higher  the  value  of  It,  the 
lower  net  heat  flow  due  to  radiation  and  convection  is  achieved.  R^t  is  the  clothing 
characteristic  that  accounts  for  water  vapor  flow  due  to  clothing  permeability.  The 
higher  the  value  of  evaporative  resistance,  the  less  evaporative  cooling  occurs;  hence,  the 
higher  the  level  of  heat  stress.  Although  It  is  associated  with  R^t,  the  relationship  is 
neither  linear  nor  fixed  for  all  clothing. 

To  complicate  matters,  dry  heat  exchange  and  evaporative  cooling  are  altered 

with  air  and  body  movement.  Consequently,  as  work  demands  increase,  it  is  possible  to 
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see  a  decrease  in  the  total  evaporative  resistance  and  IT.  Both  terms  have  static  and 
dynamic  values  associated  with  its  use.  That  is  Ii,stat  is  associated  with  the  total  insulation 
of  clothing  absent  of  movement,  and  Ij.dyn  is  associated  with  the  total  insulation  of 
clothing  with  movement.  The  same  associations  are  true  for  Re/r-  Thus,  it  is  important  to 
understand  the  thermal  resistance  properties  of  ensembles  and  how  environment  and 
activity  level  alters  them. 


Components  of  Insulation  and  Evaporative  Resistance 

In  1955  Burton  and  Edholm  introduced  the  new  unit  for  clothing  insulation  -  the 

clo.  One  clo  of  insulation  was  intended  to  be  equivalent  to  thermal  insulation  required  to 

keep  a  sedentary  person  comfortable  with  normal  indoor  clothing  at  normal  indoor 

climatic  conditions  (21°C).  The  purpose  of  using  the  unit  clo  was  to  remove  the 

awkward  physical  unit  of  m2  °C/W,  so  one  clo  equals  0. 1 55  m2  °C/W  [5],  Goldman 

points  out  the  advantage  of  using  the  clo  as  the  unit  of  insulation  is  that  it  can  be 

expressed  as  heat  loss  that  will  occur  for  the  average  adult  male  who  has  1 .8  m2  of 

surface  area,  using  a  simple  relationship  that  such  an  individual  will  lose  10  kcal/hr  of 

heat  by  radiation  and  convection  for  every  degree  (°C)  difference  between  the  average 

skin  temperature  and  the  air  temperature  with  1  clo  unit  of  insulation  [6].  Therefore  5 

# 

kcal/hr  will  be  lost  with  2  clo  units  of  insulation. 
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Total  clothing  insulation  (It)  is  the  combined  insulation  provided  by  clothing  and 
the  surrounding  layer  of  air.  Parsons  [7]  describes  this  relationship  mathematically  as: 


It  -  Id  +  la 


(2) 


Intrinsic  clothing  insulation  (Ici)  is  a  characteristic  of  the  clothing  itself  and  not 
the  external  environment  or  the  body  condition.  Ici  represents  the  resistance  to  heat 
transfer  between  the  clothing  surface  and  the  skin.  Typical  units  are  °C  m  /W.  Ici  values 
and  clo  units  are  still  used  in  several  thermal  comfort  and  clothing  standards  and 
information  on  determining  Ici  from  measured  values  of  It  is  described  in  ISO  Standard 
9920  [7,  8]. 

k 

Ia  describes  the  thermal  resistance  or  insulation  provided  by  the  air  between  the 
skin  and  garment.  The  properties  of  this  layer  are  important  to  heat  exchange  and  can  be 
affected  by  the  external  environmental  conditions. 

For  an  individual  wearing  an  ensemble,  the  surface  area  of  the  individual  is 
increased  by  an  amount  related  to  the  thickness  of  the  clothing  layer.  This  new  surface 
area  is  difficult  to  determine,  but  is  important  for  other  relationships  with  heat  transfer. 

A  clothing  adjustment  factor  (fct)  is  used  to  account  for  this  new  surface  area.  The  term 
fci  is  the  ratio  of  the  clothed  surface  area  of  the  body  to  the  nude  surface  area  of  the  body. 
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The  following  equation  is  an  approximation  for  fci  given  by  McCullough  and  Jones 
(1984)  [7], 

fci  =1.0  +  0.31  Ic,(clo)  (3) 

To  determine  the  intrinsic  clothing  insulation,  ITjStat  is  measured  using  a  clothed 
manikin  or  hot  plate  as  described  in  the  following  section.  Ia  is  measured  in  a  similar 
fashion,  but  without  the  fabric  sample  or  clothing.  Then, 


Id  I 'n  slat  ■  f, /fci  (4) 

A  more  convenient  term  for  measurement  is  effective  clothing  insulation  (Icie), 
which  is  an  approximation  for  Id  for  the  test  conditions.  This  is  described  by  the 
following  equation: 


Ide  -  It  -  la  (5) 

This  same  principle  can  be  applied  for  the  total  evaporative  resistance  of  clothing 
(Re/r)  by  dividing  it  into  two  components.  The  evaporative  resistance  due  to  the  clothing 
itself  (Rd)  and  that  due  to  the  air  layer  (Ra)  near  the  clothing  or  exposed  skin. 


Re.t  -  Rd  +  Ra 
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(6) 


Values  for  Re>T  and  Ra  can  be  determined  empirically  from  variations  of  the 
standard  tests  for  clothing  insulation  using  sweating  hot  plates  or  sweating  manikins  (see 
next  section).  In  this  way,  Rci  can  be  estimated  from  the  following  equation. 


RcI  =  Re,T-Ra/fcl  (7) 

Again,  a  convenient  approximation  is 

Rcfe  —  Re,T  "  Ra  (^0 


Laboratory  Test  Methods 

There  are  three  different  methods  for  determining  the  thermal  properties  of  a 
garment.  The  first  method  involves  the  use  of  a  heated  plate;  the  second  involves  a 
heated  copper  manikin;  and  the  third  method  involves  the  use  of  human  participants. 


Hot  Plate  Method 

While  manikins  and  hotplates  use  similar  basic  principles  to  determine  heat  loss 
and  insulation  values,  they  typically  have  different  end  goals.  A  hotplate  is  designed  to 
provide  accurate  one-dimensional  heat  and  moisture  flow  through  a  fabric  sample  to 
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determine  thermal  and  water  vapor  resistance.  The  goal  is  to  accurately  evaluate  the 
material  properties  for  the  test  environment  only  [9]. 

The  American  Society  for  Testing  and  Materials  (ASTM)  developed  a  standard 
method  for  using  a  sweating  hot  plate  in  method  F  1868-02  [10].  This  test  method  covers 
the  measurement  of  the  thermal  resistance  and  the  evaporative  resistance  under  steady- 
state  conditions  of  fabrics,  films,  coatings,  foams,  and  leathers,  including  multi-layer 
assemblies,  for  use  in  clothing  systems.  There  are  several  relevant  measures  from 
sweating  hot  plate  tests.  The  most  basic  measure  is  the  operating  heat  flux  required  to 
maintain  a  constant  skin  temperature.  In  dry  tests,  it  represents  the  conductive/ 
convective/radiative  heat  transfer.  In  sweating  tests,  it  also  includes  evaporative  heat 
losses.  This  sweating  test  is  the  most  common  method  used. 


Copper  Manikin 

A  life-sized  heated  copper  manikin  can  be  used  in  the  evaluation  of  the  heat 
transfer  potential  of  clothing  garments.  Similar  to  the  hot  plate  method,  the  manikin  is 
electrically  heated  so  that  the  skin  temperature  is  similar  to  that  of  people.  Depending  on 
number  of  surface  segments  the  resolution  can  be  adjusted  to  be  sufficiently  high  to 
complete  the  measurement  task.  Some  manikins  in  use  today  have  1  zone  while  others 
have  more  than  100  individually  regulated  segments.  By  summing  up  the  area  weighted 
heat  loss  values  from  the  manikin,  a  total  value  for  whole  body  heat  loss  is  determined. 
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Some  performance  features  of  the  most  commonly  used  thermal  manikins  are: 
simulation  of  human  body  heat  exchange,  measurement  of  3 -dimensional  heat  exchange, 
integration  of  dry  heat  losses,  measurement  of  clothing  thermal  insulation,  product 
development,  and  providing  values  for  prediction  models. 

Values  obtained  with  different  manikins  in  different  laboratories  should  be 
comparable  and  similar  within  defined  limits  for  the  same  test  conditions.  The  conditions 
and  requirements  for  comparable  measurements  with  different  manikins  and  methods  are 
defined  in  standards.  The  American  Society  for  Testing  and  Materials  standardized  this 
procedure  in  ASTM  method  F  1291-99,  and  the  International  Organization  for 
Standardization  standardized  the  procedures  in  ISO  9920  [11, 12].  ASTM  method  F 
1291  and  ISO  9920  are  both  in  the  process  of  being  updated  to  reflect  the  changes  with 
the  new  sweating  manikin.  Similar  to  the  hot  plate  procedures,  these  test  methods  cover 
the  measurement  of  the  thermal  resistance  and  the  evaporative  resistance  under  steady- 
state  conditions.  With  over  a  hundred  different  manikins  being  built  and  used  around  the 
world,  it  is  difficult  for  any  standard  to  encompass  procedures  for  all  types  of  manikins. 

Thermal  manikins  have  evolved  from  its  first  model  in  1941  for  testing  military 
clothing,  and  can  be  grouped  into  three  categories.  First  are  static  (non-moving)  and  non¬ 
perspiring  units,  second  are  movable  (walkable),  but  non-perspiring  ones  such  as  the 

copper  manikin  ‘Charlie’  in  Germany  used  by  Mecheels  and  Umbach  in  1977,  and  third 
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are  sweating  manikins  [13].  To  simulate  sweating  on  non-perspiring  manikins,  many 
researchers  used  highly  absorbent  fabrics  on  the  manikin  (under  the  tested  garment)  and 
supplied  water  to  the  “underwear”  by  sprinkling  or  water  pipes.  The  third  generation 
manikins  simulate  true  perspiration  and  body  motion,  but  again  not  all  sweating  manikins 
are  dynamic. 

Recent  innovation  at  the  Institute  of  Textiles  and  Clothing,  Hong  Kong,  produced 
Walter  the  "sweating"  manikin.  Walter  is  made  up  of  water,  mechatronics  and  breathable 
fabric,  allowing  realistic  simulation  of  human  thermal  physiology  under  various 
environments.  Walter  has  waterproof  but  moisture-permeable  fabric  skin,  which  can  be 
unzipped,  removed  and  interchanged  with  different  skins,  simulating  different  rates  and 
patterns  of  perspiration.  In  addition,  Walter  is  a  dynamic  manikin  that  allows  researchers 
to  simulate  the  process  of  walking.  This  new  technology  surprisingly  has  an  affordable 
price  tag  associated  with  it  as  it  may  be  90%  less  than  traditional  copper  and  plastic 
manikins  [14]. 

Human  Tests  for  Clothing  Thermal  Characteristics 

Though  useful,  measurements  of  clothing  thermal  properties  made  on  hot  plates 
or  manikins  do  not  represent  the  properties  of  clothing  during  wear.  The  movements  of 
the  worker  increase  the  convective  air  flow  both  between  layers  and  at  the  clothing 
surface,  modifying  both  insulation  and  vapor  permeability.  Although  recent 
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developments  with  dynamic  sweating  manikins  provide  accurate  data  on  clothing 
insulation  and  evaporative  resistance,  there  still  exists  a  difference  between  manikin  and 
human  wear  tests. 

Human  wear  testing  to  determine  evaporative  resistance  and  clothing  insulation 
values  are  based  on  determining  the  “prescriptive  zone”  as  described  by  Lind  [15].  The 
upper  limit  prescriptive  zone  was  defined  as  the  point  where  heat  loss  and  heat  gain  were 
equal  (Emgx —  Ereq). 

Using  the  premise  of  the  prescriptive  zone,  Belding  and  Kamon  proposed  a 

method  for  determining  ambient  vapor  pressures  at  36°C  for  a  variety  of  exercise 

intensities  and  air  movements  [16].  Their  method  used  a  time-intensive  protocol  to 

determine  critical  environmental  conditions  for  evaporative  heat  loss.  Kamon  and 

Avellini  used  the  same  approach  as  Belding  and  Kamon  [17].  In  their  experiments,  the 

participants  were  subjected  to  a  range  of  ambient  temperatures  between  36  and  52°C  with 

the  water  vapor  pressure  progressively  increased  at  each  ambient  temperature.  It  was 

expected  that  on  the  basis  of  the  body  core  temperature  inflection,  a  line  for  the  safe  limit 

for  the  psychrometric  chart  would  be  empirically  identified  particularly  for  the  higher 

temperatures,  where  individual  sweating  rate  capacity  was  believed  to  be  the  limiting 

factor.  Holmer  and  Elnas  developed  a  method  to  determine  both  the  evaporative  and 

sensible  heat  loss  occurring  simultaneously  [18].  However,  this  method  was  difficult 

because  direct  measurement  of  the  water  vapor  pressure  gradient  between  the  skin  and 
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ambient  air  was  required.  Kenney  et  al.  simplified  this  procedure  by  minimizing  the 
number  and  duration  of  tests  necessary  to  determine  these  limits  [19].  However,  these 
more  time-efficient  protocols  require  the  ability  to  systematically  change  ambient 
temperature  or  water  vapor  pressure. 

In  the  Kenney  method,  the  metabolic  rate  target  was  30%  of  maximal  aerobic 
capacity  (MAC).  Generally,  a  person  can  work  at  1/3  their  MAC  for  an  8-hour  work  day 
[20]. 


The  testing  chamber  was  controlled  in  that  the  dry  bulb  temperature  (Tab)  and  wet 
bulb  temperature  (TWb)  could  be  closely  manipulated.  Air  velocity  was  0.5  m/s  or  less. 
The  participant’s  heart  rate  (HR),  rectal  temperature  (Tre)  and  mean  skin  temperature 
(TSk)  were  monitored. 

Each  participant  partook  in  two  trials  wearing  the  garment  to  be  tested.  In  one  of 
the  trials,  Tab  was  held  constant  and  after  30  minutes  for  stabilization,  the  ambient  water 
vapor  pressure  (Pa)  was  increased  in  increments  of  0.13  kPa  every  5  minutes.  In  the  other 
trial,  the  Pa  was  held  constant  at  a  low  humidity  while  the  Tab  was  increased  in  1°C 
increments  every  five  minutes. 

The  Tre  for  each  participant  was  plotted  and  the  point  of  inflection  where  Tre 

sharply  rose  was  noted.  This  inflection  point  represented  the  inability  of  the  body  to 
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dissipate  the  heat  load  and  thereafter  heat  was  stored  by  the  body.  Data  from  a 
participant  was  plotted,  as  shown  in  Figure  1,  to  illustrate  the  typical  time  course  for  an 
inflection  point  protocol.  The  important  points  of  this  figure  is  the  starting  Tre  with  a 
steady  rise  until  a  steady  state  is  achieved,  and  then  the  inflection  point  followed  by  steep 
rise  in  Tre.  Using  this  method,  the  critical  temperature  (Tcrit)  at  a  given  Pa  and  the  critical 
water  vapor  pressure  (Pent)  at  a  given  Tab  can  be  determined. 


Time  (min) 


Figure  1 .  Typical  Time  Course  for  Tre  During  an  Inflection  Point  Protocol. 


At  the  inflection  point,  the  required  rate  of  evaporative  cooling  (Ereq)  was  equal  to 
the  maximum  rate  of  evaporation  (Emax).  In  other  words,  the  rate  of  heat  storage  was  zero 
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since  the  evaporative  cooling  was  equal  to  the  net  heat  gain  from  metabolism  plus  the  dry 
heat  exchange. 

Emax  =  (Mnet)  +  (R+C)  (9) 

Another  condition  that  existed  at  the  inflection  point  was  that  evaporative  cooling 
was  at  its  maximum  value.  Here  evaporative  cooling  (Emax)  was  equal  to  the  difference 
in  water  vapor  pressure  between  the  skin  (PSk)  and  the  ambient  environment  (Pa)  divided 
by  the  total  evaporative  resistance  (Rej).  This  is  shown  in  the  following  equation: 

Emax  —  (Psk"Pa)/Re.T  (10) 


The  relationship  of  Re/r  with  respect  to  AP  (PSk  -  Pa)  and  Emax  is  based  on  the 
assumption  that  R^t  remains  the  same  as  AP  and  Emax  change.  However,  there  are  some 
researchers  that  question  this  principle.  In  his  paper,  Bernard  found  that  a  warm  humid 
environment  resulted  with  a  lower  R^t  [30].  Theoretically,  the  maximum  sweat  rate  is 
proportional  to  the  difference  in  the  saturated  partial  pressure  of  water  at  the  skin  minus 
the  partial  pressure  in  the  air,  and  the  evaporative  resistance  of  clothing  worn  has  no 
effect.  The  question  about  R^t  being  a  constant,  as  accepted  in  Equation  10,  regardless  of 
environment  has  not  been  evaluated. 
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Without  a  direct  source  of  radiant  heat,  the  rate  of  dry  heat  exchange  (R+C)  was 
taken  as  the  difference  in  Tab  and  TSk  divided  by  the  clothing  insulation  (It).  This  is 
presented  in  the  following  equation: 

(R+C)  =  (Tdb-Tsk)/IT  (11) 

By  substituting  Equations  10  and  1 1  Into  Equation  9,  the  following  equation  results. 

(Psk"Pa)/Re,T  —  (Mnet)“^(Tab-Tsk)/lT  (12) 

When  the  measured  values  and  environmental  conditions  were  placed  in  Equation 
12  for  each  of  two  inflection  points,  there  were  two  equations  with  two  unknowns.  This 
allowed  for  the  calculation  of  lx  and  Re/r. 

At  each  inflection  point,  heat  gain  equals  heat  loss.  This  is  represented 
mathematically  using  the  following  equation: 

Mnet  +  (R  +  C)  =  E,  (13) 

where  Mnet  is  the  net  metabolic  heat  production  (M)  corrected  for  external  work  (W)  and 
respiratory  exchanges  due  to  convection  (Cres)  and  evaporation  (Eres).  Metabolic  rate  (M) 
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in  W/m2  was  estimated  from  oxygen  consumption  in  liters  per  minute  and  the  respiratory 
ratio  (R)  using  the  following  equation  [21]: 

M  =  352(0.23R  +  0.77)  V02/Ad  (14) 

The  Dubois  surface  area  (Ad)  was  calculated  for  each  subject  using  the  following 
equation  [22]: 


AD  =  0.202  W° 425  H0725,  (15) 


where  W  was  the  weight  of  the  body  (kg)  and  H  was  the  height  of  the  body  (m). 


The  external  work  (W)  was  calculated  (W/m2)  using  the  following  equation: 


W  —  0.163  nib  Vw  fg/Ao,  (16) 

where  mb  was  body  mass  in  kg,  Vw  was  walking  velocity  in  m/min,  fg  was  the  fractional 
grade  of  the  treadmill,  and  Ad  was  the  Dubois  surface  area. 

Respiratory  exchanges,  latent  respiration  heat  loss  (Eres)  and  dry  respiration  heat 
loss  (Cres),  were  calculated  as  follows: 
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Cres  =  0.0012  M  (34-Tdb) 


(17) 


and 


Eres  =  0.0173  M  (5.87-Pdp) 


(18) 


^5 

The  net  metabolic  rate  (Mnet)  from  Equation  13  can  be  calculated  in  W/m  using 
the  following  equation: 


Mnet  =  (M-W)  +  Cres- Eres  (19) 

Kenney  et  al.  recognized  that  there  may  be  some  heat  storage  represented  by  a 
gradual  change  in  Tre  [21].  To  account  for  this,  the  rate  of  change  in  heat  storage  can  be 
estimated  knowing  the  specific  heat  of  the  body  (0.97  W  h/°C  kg),  body  weight  (BW), 
and  the  rate  of  change  of  body  temperature  (ATre  /  At)  before  the  inflection  point  was 
reached  [23].  That  is, 

S  =  0.97  BW  ATre  /  Ad  At  (20) 
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Thermal  Resistance  Values  for  Various  Work  Ensembles 

By  using  a  hot  plate,  copper  manikin  or  human  subjects,  thermal  resistance  values 
for  garment  ensembles  can  be  quantified  in  terms  of  total  insulation  (IT)  and  total 
evaporative  resistance  (R^t).  Experimentally  determined  values  for  select  work 
ensembles  are  presented  in  Table  1  for  comparison  purposes. 


The  It  and  Re  T  values  reported  in  Table  1  vary  within  each  garment  between  the 
method  used.  Havenith  et  al.  found  that  heated  manikin  results  for  standing/no  wind 
appear  to  be  on  the  average  0.15  clo  (0.023  °C  m2/W)  higher  than  human  subjects  results 
[25]. 


Table  1,  Comparison  of  Experimentally  Determined  Thermal  Resistance  Values. 


Ensemble  Description 

27^ 

O 

o 

Re/r 

(kPa  m2/W) 

Reference 

Method 

Tyvek®  Coverall 

0.070 

0.020 

[24] 

hot  plate 

Tyvek®  Coverall 

0.171 

0.033 

[24] 

manikin 

Tyvek®  Coverall 

0.086 

0.017 

[23] 

human 

participant 

Gore-Tex®  Outer-wear 

0.054 

0.009 

[24] 

hot  plate 

Gore-Tex®  Outer-wear 

0.210 

0.032 

[24] 

manikin 

Gore-Tex®  Outer-wear 

0.130 

0.028 

[21] 

human 

participant 

Cotton,  Single  Knit 

0.079 

0.009 

[24] 

hot  plate 

Men’s  Summer  Casual 
(short  sleeve) 

0.201 

0.029 

[24] 

manikin 

Military  Fatigues 

0.090 

0.016 

[21] 

human 

participant 
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The  data  in  Table  1  shows  a  larger  gap  (0.52  -  0.72  clo)  between  human 
participants  and  manikin  data.  Nishi  et  al.  [26]  and  Vogt  et  al.  [27]  support  Havenith’s 
findings,  but  Nielsen  at  al.  [28]  and  Olesen  et  al.  [29]  found  the  human  participant  values 
were  0.22  clo  lower  than  manikin  data. 

Havenith  et  al.,  found  that  the  permeability  index  (I’m)  changed  with  wind  and 
movement  [30],  In  their  study,  the  permeability  increased  three  fold  with  permeable 
clothing  and  six  fold  with  impermeable  clothing.  Additionally,  they  found  that  the  total 
insulation  was  reduced  by  32%  and  that  walking  at  slower  rates  yielded  smaller  gains. 
Breckenridge  and  Goldman  [31]  reported  similar  findings  with  an  increase  in  Im  by  123% 
and  a  decrease  in  It  by  28%. 

A  few  years  later,  Parsons  et  al.  [32],  Hoimer  et  al.  [5],  and  Havenith  at  al.  [33]  all 
find  that  It, stat  needs  to  be  adjusted  for  wind  and  walking.  Their  findings  were  adapted  in 
ISO  7933.  Havenith  et  al.  reports  that  a  78%  reduction  of  Rejwith  a  50%  reduction  in  It 
can  be  seen  [33]  which  is  similar  to  the  differences  seen  in  Table  1 .  I^dyn  is  converted  by 
multiplying  Ir,statby  a  correction  factor  (CFci)  as  shown  in  Equations  21  and  22: 


I'l'.dvn  CFci  X  It, stat 


(21) 


and 
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(22) 


pr?  _  (o.043-0.39SVar+0.66V^r-0.37SWaIksp+0.094Walksp2) 

K-'rcl  ~  e 

where  Var  is  the  velocity  of  the  air  and  Walksp  is  the  walking  speed.  The  speed  is 
calculated  based  on  the  metabolic  demand  (M  in  W/m2)  and  is  shown  in  the  equation 
below: 


Walksp  =  0.0052  (M  -  58)  (23) 

Another  obvious  relationship  in  Table  1  is  that  the  hot  plate  values  are  less  than 
the  manikin  values  for  both  It  and  Rej.  There  are  a  few  possible  reasons  for  these 
differences. 

First,  the  values  may  be  lower  than  manikin  data  because  of  the  fit  or  drape  of  the 
clothing.  Hot  plates  tend  to  tested  with  a  tight  fit  where  the  manikins  and  humans  use  a 
looser  fit.  Second,  the  wetting  of  the  material  is  likely  to  alter  the  values.  Manikin 
values  were  measured  on  dry  manikins  to  determine  the  dry  heat  exchange  while  the  hot 
plates  were  wet.  The  clothing  on  human  participants  was  also  wet  from  sweating. 
Although  the  hot  plate  values  are  not  the  same  as  that  of  the  human  participants,  they  are 
close  in  two  of  the  garments  tested.  Wetting  of  the  clothing  alters  the  obtained  value  by 
attenuating  the  resistance  [34].  Third,  it  is  difficult  to  simulate  the  movements  of 
exercising  humans.  The  presence  of  body  motion  aids  in  the  circulation  of  air  through 
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the  clothing  and  therefore  also  reduces  the  resistance.  The  effects  of  air  and  body 
movement  on  It  and  Rej  have  been  well  documented  [5, 21, 23, 21,  and  28].  Although 
there  is  agreement  of  this  needed  adjustment,  manikin  data  is  not  always  adjusted  for  air 
and  movement. 

The  methodology  for  using  human  participants  in  conducting  heat  stress  studies 
proposed  by  Kenney  et  al  [19]  has  been  used  by  Bernard  and  Matheen  [9],  Barker  et  al. 
[23],  Kenney  and  Zeman  [35],  and  Malcolm  et  al.  [36]  as  well  as  other  researchers. 


Hypothesis 

A  reasonable  evaluation  of  selected  protective  clothing  garments  would  be  a 
determination  of  their  heat  exchange  characteristics.  The  primary  purpose  of  this  paper 
was  to  explore  the  methodology  for  being  able  to  distinguish  between  garments  based  on 
the  total  evaporative  resistance  properties  within  different  environments  and  work 
demands.  The  secondary  purpose  is  to  challenge  the  relationship  of  R^t  with  respect  to 
changes  in  AP  and  Emax.  The  default  assumption  is  that  R^t  remains  the  same  as  AP  and 
Emax  change. 

There  are  three  null  hypothesis  to  be  tested:  (1)  there  are  no  differences  between 

mean  Rej  values  among  ensembles,  (2)  there  are  no  differences  between  mean  Revalues 

among  environments  and  metabolic  rates/demands,  and  (3)  there  are  no  differences 
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between  mean  Re,T  values  while  AP  changes  within  environments  and  metabolic 
demands. 
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METHODS 


The  primary  purpose  of  this  paper  was  to  explore  the  methodology  for 
distinguishing  between  garments  based  on  the  total  evaporative  resistance  properties 
within  different  environments  and  work  demands.  The  secondary  purpose  is  to  challenge 
the  relationship  of  Re,i  with  respect  to  changes  in  AP  and  Emax.  Experimental  trials  were 
conducted  to  determine  the  evaporative  resistance  for  five  clothing  ensembles.  The 
protocols  included  a  fixed  metabolic  demand  under  three  different  relative  humidity 
levels  for  Phase  1,  and  for  three  metabolic  demands  with  a  fixed  relative  humidity  level 
for  Phase  2.  The  key  to  these  studies  was  being  able  to  distinguish  the  point  of  transition 
from  compensable  heat  stress  to  uncompensable  heat  stress  (Ereq  =  Emax). 


Participants 

Fourteen  adults  (nine  men  and  five  women)  participated  in  experimental  trials  for 
Phase  1  and  fifteen  adults  (eleven  men  and  four  women)  participated  in  Phase  2.  Their 
physical  characteristics  are  provided  in  Appendix  A  and  the  average  and  standard 
deviation  of  their  physical  characteristics  by  gender  are  provided  in  Table  2,  Following 
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the  local  IRB  procedures,  a  written  informed  consent  was  obtained  and  all  subjects  were 
qualified  by  a  physician. 

Prior  to  beginning  the  experimental  trials,  participants  underwent  a  5 -day 
acclimatization  period.  Acclimatization  involved  walking  on  a  treadmill  at  a  metabolic 
rate  of  approximately  160  W/m2  in  a  climatic  chamber  at  50°C  and  20%  relative 
humidity  (rh).  During  acclimation  participants  wore  shorts  (and  sports  bra),  socks  and 
shoes. 


Table  2.  Summary  of  Participant  Characteristics. 


Protocol 

Gender 

Num 

Age 

(yrs) 

Height 

(cm) 

Weight 

(kg) 

Surface 
Area  (m2) 

Phase  1 

Men 

9 

29.2  ±6.8 

183  ±6.0 

97.2  ±  18.5 

2.18  ±0.20 

Humidity 

Women 

5 

31.8  ±  9.1 

161  ±7.0 

63.5  ±  17.2 

1.66  ±0.23 

All 

14 

30.1  ±7.5 

175  ±  12.0 

85.2  ±24.1 

2.00  ±  0.33 

Pjitiep  9 

Men 

11 

28.0  ±9.5 

176  ±  11.2 

81.9  ±11.7 

1.98  ±0.18 

Metabolic 

Women 

4 

23.0  ±4.7 

165  ±6.3 

64.2  ±  18.0 

1.70  ±0.22 

All 

15 

26.7  ±8.6 

173  ±  11.1 

77.2  ±  15.3 

1.91  ±0.22 

Men 

20 

28.6  ±8.2 

180  ±8.6 

89.55  ±  15.1 

2.08  ±0.19 

Both 

Women 

9 

27.4  ±6.9 

163  ±6.7 

63.85  ±17.6 

1.68  ±0.23 

All 

29 

28.0  ±7.5 

171  ±7.6 

76.7  ±  16.35 

1.88  ±0.21 

Clothing  Ensembles 

Five  different  clothing  ensembles  were  evaluated  in  each  Phase  with  only  one 
ensemble  being  changed  for  Phase  2.  The  ensembles  included:  Ensemble  A  —  work 
clothes  (4  oz/yd2  cotton  shirt  and  8  oz/yd2  cotton  pants);  Ensemble  B  —  cotton  coveralls 
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(9-10  oz/yd2 )  and  three  limited-use  protective  clothing  ensemble:  Ensemble  C  --  particle- 
barrier  ensembles  (Tyvek®  1424  for  Phase  1  and  Tyvek  1427  for  Phase  2),  Ensemble  D  - 
-  water-barrier,  vapor-permeable  ensembles  (NexGen®  LS  417),  and  Ensemble  E  — 
vapor-barrier  ensembles  (Tychem®  QC).  The  limited-use  coveralls  had  a  zippered 
closure  in  the  front  and  elastic  cuffs  at  the  arms  and  legs. 

All  ensembles  were  worn  without  a  hood  and  a  cotton  tee-shirt  and/or  sports  bra 
and  shorts  were  worn  under  all  clothing  ensembles. 


Protocols 

Three  experimental  protocols  were  followed  each  Phase.  The  design  for  Phase  1 
had  three  environments  with  a  fixed  metabolic  rate.  Treadmill  speed  and  grade  were  set 
to  elicit  a  metabolic  rate  of  about  160  W/m2.  The  first  protocol  (R7)  was  a  warm/humid 
environment  designed  to  reduce  Emax  by  limiting  evaporation.  The  second  protocol  (R2) 
was  a  hot/dry  environment  designed  to  increase  Ereq  by  increasing  radiant  and  convective 
(R+C)  heat  gains.  The  third  protocol  (R5)  was  a  moderate  environment  designed  to 
increase  R+C  while  decreasing  Emax. 

In  the  R7  protocol,  the  dry  bulb  temperature  (Tdb)  was  set  at  30°C  and  relative 
humidity  (rh)  at  70%.  Once  the  participant  reached  thermal  equilibrium  (no  change  in  Tre 
and  heart  rate  for  at  least  15  minutes),  Tab  was  increased  0.7°C  every  5  minutes.  In  the 
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R2  protocol,  Tdb  was  set  at  40°C  with  rh  at  20%.  When  participants  reached  thermal 
equilibrium,  Tdb  was  increased  1°C  every  5  minutes.  For  the  R5  protocol,  Tdb  was  set  at 
34°C  with  50%  rh.  Upon  reaching  thermal  equilibrium,  Tdb  was  increased  0.8°C  every  5 
minutes. 

In  Phase  2,  the  study  design  called  for  three  metabolic  demands:  light  demand, 
metabolic  rate  of  80  W/m2  (Ml);  moderate  demand,  metabolic  rate  of  160  W/m2  (M2); 
and  heavy  demand,  with  a  metabolic  rate  of  240  W/m2  (M3).  Actual  metabolic  rates 
were  calculated  using  oxygen  consumption  based  on  open  circuit  indirect  calorimetry  and 
body  surface  area. 

The  environment  was  set  at  a  50%  relative  humidity  (rh).  The  starting 
temperature  for  the  trials  was  set  at  34°C,  but  varied  based  on  the  ensemble  being  worn 
and  individual.  When  participants  reached  thermal  equilibrium,  Tdb  was  increased  1°C 
every  5  minutes. 

Trials 

The  trials  were  conducted  in  a  Model  7010  climatic  chamber  designed  by  Forma 
Scientific.  The  chamber  was  2.4  m  wide,  3.0  m  deep,  and  2.2  m  high  (8.0  x  10.0  x  7.3 
ft).  The  range  of  humidity  was  10-90%  and  the  temperature  range  was  4-60°C  (40- 
140°F).  Temperature  and  humidity  were  controlled  according  to  protocol  and  air  speed 
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was  0.5  m/s.  The  work  demand  consisted  of  walking  on  a  motorized  treadmill  at  a  speed 
and  grade  set  to  elicit  the  desired  metabolic  rate  (80, 160,  or  240  W/m2). 

Heart  rate  was  monitored  using  a  Polar  heart  rate  monitor.  Core  temperature  was 
measured  with  a  flexible  YSI  thermistor  (401AC)  inserted  10  cm  beyond  the  anal 
sphincter  muscle.  The  thermistor  was  calibrated  prior  to  each  trial  using  a  hot  water  bath. 
Skin  temperatures  were  measured  with  an  YSI  surface  thermistor  (409AC)  taped  to  the 
skin  at  four  points  (left  chest,  right  upper  arm,  right  thigh,  and  left  calf).  Average  skin 
temperature  was  determined  using  a  modified  Ramanathan  Technique  as  shown  in  the 
following  equation  [7]: 


Tsk  0.3  Tchest  "t"  0.3  Tantt  0.2  Tthigh  0.2  Tcalf 


(25) 


Assessment  of  oxygen  consumption  was  used  to  establish  metabolic  rate. 
Participants  breathed  through  a  two-way  valve  connected  to  flexible  tubing  that  was 
connected  to  a  collection  bag  (Douglas  bag).  Expired  gases  were  collected  every  30 
minutes  during  the  experiments  for  2.5  minutes.  The  volume  of  expired  air  was 
measured  using  a  dry  gas  meter,  A  small  aliquot  was  removed  from  the  Douglas  bag  and 
drawn  through  a  drying  agent  (DriRite)  into  a  Beckman  Model  E2  Oxygen  Analyzer  to 
determine  oxygen  content.  Oxygen  consumption  (VO2)  was  calculated  according  to 
Equation  25. 
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vo2  =  vE  ao2  cf 


(26) 


Where,  VE  was  the  expired  air  flow  rate  in  liters  per  minute,  A02was  the 
difference  in  the  fraction  of  oxygen  between  the  inspired  and  expired  air,  and  CF  was  a 
correction  factor  to  convert  the  volume  to  standard  temperature  and  pressure  dry  (STPD) 
[35]. 

During  trials,  participants  were  allowed  to  drink  water  or  a  commercial  fluid 
replacement  beverage  at  will. 

Core  temperature,  heart  rate  and  ambient  conditions  (dry  bulb,  psychrometric  wet 
bulb  and  globe  temperatures)  were  monitored  continuously  and  recorded  every  5 
minutes.  Trials  lasted  approximately  120  minutes  unless  one  of  the  following  criteria 
was  met:  (1)  a  clear  rise  in  Tre  associated  with  a  loss  of  thermal  equilibrium,  (2)  Tre 
exceeded  39  °C,  (3)  a  sustained  heart  rate  greater  than  85%  of  the  age-predicted 
maximum  heart  rate,  or  (4)  participant  wished  to  stop. 

The  order  of  the  ensemble-environment  conditions  was  randomized.  Any  trial 
that  had  to  be  repeated  was  repeated  at  the  end.  An  experimental  trial  data  dictionary  is 
presented  with  the  data  for  each  phase  in  the  appendices.  Phase  1  data  are  provided  in 
Appendix  B  and  Phase  2  data  are  in  Appendix  C. 
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Critical  Conditions 


By  evaluating  the  point  at  which  a  clear  rise  in  Trej  associated  with  a  loss  of 
thermal  equilibrium,  the  critical  condition  (Ereq  =  Emax)  can  be  determined  by  using  the 
data  point  preceding  this  rise.  At  the  point  of  critical  conditions  R^t  can  be  calculated 
using  Equation  12.  In  this  equation  there  are  two  unknowns,  Re^dyn  and  I-r.dyn-  The  IT>stat 
values  were  calculated  from  measured  insulation  values  (clo)  according  to  ASTM  F 
1291,  Standard  Test  Method  for  Measuring  the  Thermal  Insulation  of  Clothing  using  a 
Heated  Manikin,  Option  #1  [12]. 

The  insulation  provided  by  clothing  (clo)  was  measured  using  an  electrically- 
heated  manikin  in  thermal  equilibrium  with  the  surrounding  environment.  The  manikin 
is  a  full  size  male  with  19  electrically  separate  segments.  The  manikin  has  knee,  hip, 
shoulder,  and  elbow  joints  that  can  be  flexible  or  locked  in  an  immobile  position  [11]. 

Measurement  and  control  of  the  heat  supply  for  each  section  is  achieved  by  using 
a  digital  process  computer.  Display  and  recording  of  the  data  is  conducted  by  a  second 
computer  which  is  serially  interfaced  with  the  process.  Temperature  readings  and  power 
input  values  for  each  segment  are  area  weighted  when  calculating  the  total  insulation 
value  [11]. 
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The  insulation  value  (clo)  was  measured  according  to  ASTM  F  1291,  Standard 
Test  Method  for  Measuring  the  Thermal  Insulation  of  Clothing  using  a  Heated  Manikin, 
Option  #1  [11].  The  chamber  had  an  ambient  air  temperature  of  20°C,  dew  point 
temperature  was  controlled  at  1°C  and  air  velocity  of  0.2  m/s,  and  the  manikin  surface 
temperature  was  set  at  33.2°C. 

To  test  each  ensemble,  the  manikin  was  dressed  in  an  ensemble  with  all  closures 
secured.  It  was  hung  from  a  metal  stand  by  a  hook  in  the  head.  The  feet  touched  the 
floor  with  the  arms  hung  at  the  sides.  Equilibrium  was  maintained  for  at  least  one  hour 
prior  to  testing.  Data  were  collected  by  computer  every  30  seconds  for  the  30  minute  test 
[11]. 


The  Ij.dyn  values  were  calculated  for  each  ensemble  by  adjusting  the  IT>stat  values 
for  wind  and  speed  as  suggested  by  Havenith  et  al  [32]. 

Additionally,  measured  trial  data  was  used  to  compute  other  variables  (AT,  AP 
and  M)  needed  to  compute  Rej.  The  metabolic  rate  was  computed  based  on  O2 
consumption  using  Equations  14,  15  and  26.  The  equations  for  differences  in 
temperatures  and  partial  pressures  are  shown  below. 

AT  =  Tab  -  Tsk  (27) 
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AP  =  Psk  -  Pi 


(28) 


Psk  =0.6105  xe{  /fc+2373«  (29) 

f(17-27xWr  +23730-o.o67x(rt)i-r  „)) 

P,=  0.6105  xe11  /(r"+237'3)J  J  (30) 


where  AP  is  the  difference  in  partial  pressure  of  water  vapor  between  the  skin  and 
ambient  air. 


Statistical  Analysis 

Statistical  analysis  included  general  descriptive  statistics  and  linear  modeling. 
The  primary  data  analysis  was  conducted  with  analysis  of  variance  (ANOVA)  and 
verified  with  the  Mixed  Procedure.  If  a  significant  difference  among  ensembles  was 
found  at  a  =  0.05,  Tukey’s  Honestly  Significantly  Different  (HSD)  was  calculated  [36]. 
If  the  difference  between  any  treatment  mean  value  was  greater  than  the  HSD,  then  the 
difference  was  determined  to  be  statistically  different. 

The  data  were  reviewed  for  outliers  defined  as  data  points  exceeding  the  mean  ± 
two  times  the  standard  deviation.  A  Sharpio-Wilkes  statistical  test  for  fit  was  performed 
on  the  ensemble  datasets  to  determine  the  best  fit  of  the  data  (normal  or  log  normal).  All 
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of  the  data  fit  well  as  being  normally  distributed.  The  data  (participant,  ensemble, 
protocol,  and  Re  t)  were  imported  into  SAS  version  8.2. 

Since  the  data  were  not  balanced,  the  data  was  analyzed  using  a  mixed  linear 
model  as  well  as  the  standard  liner  model  (GLM)  for  comparison.  The  mixed  procedure 
fits  a  variety  of  mixed  linear  models  to  data  and  enables  these  fitted  models  to  make 
statistical  inferences  about  the  data.  A  mixed  linear  model  is  a  generalization  of  the 
standard  linear  model  used  in  the  GLM  procedure,  the  generalization  being  that  the  data 
are  permitted  to  exhibit  correlation  and  non-constant  variability.  The  mixed  linear  model, 
therefore,  provides  the  flexibility  of  modeling  not  only  the  means  of  the  data  (as  in  the 
standard  linear  model)  but  their  variances  and  covariances  as  well  [37]. 

The  primary  assumptions  underlying  the  analyses  performed  by  SAS  model  Proc 
Mixed  (PM)  are  as  follows:  the  data  are  normally  distributed  (Gaussian),  the  means  of 
the  data  are  linear  in  terms  of  a  certain  set  of  parameters,  the  variances  and  covariances 
of  the  data  are  in  terms  of  a  different  set  of  parameters,  and  they  exhibit  a  structure 
matching  one  of  those  available  in  PM  [37]. 

The  fixed-effects  parameters  are  associated  with  known  explanatory  variables,  as  in 

the  standard  linear  model.  These  variables  can  be  either  qualitative  (as  in  the  traditional 

analysis  of  variance)  or  quantitative  (as  in  standard  linear  regression).  However,  the 

covariance  parameters  distinguish  the  mixed  linear  model  from  the  standard  linear  model. 

37 


The  need  for  covariance  parameters  arises  quite  frequently  in  applications.  The  most 
typical  scenarios  include:  (1)  the  experimental  units  on  which  the  data  are  measured  can 
be  grouped  into  clusters,  and  the  data  from  a  common  cluster  are  correlated,  and  (2) 
repeated  measurements  are  taken  on  the  same  experimental  unit,  and  these  repeated 
measurements  are  correlated  or  exhibit  variability  that  changes. 

PM  provides  a  variety  of  covariance  structures  to  handle  the  previous  two  scenarios. 
The  most  common  of  these  structures  arises  from  the  use  of  random-effects  parameters, 
which  are  additional  unknown  random  variables  assumed  to  impact  the  variability  of  the 
data.  The  variances  of  the  random-effects  parameters,  commonly  known  as  variance 
components,  become  the  covariance  parameters  for  this  particular  structure.  Traditional 
mixed  linear  models  contain  both  fixed-  and  random-effects  parameters,  and,  in  fact,  it  is 
the  combination  of  these  two  types  of  effects  that  led  to  the  name  mixed  model.  Proc 
Mixed  fits  not  only  these  traditional  variance  component  models  but  numerous  other 
covariance  structures  as  well. 

PM  computes  several  different  statistics  suitable  for  generating  hypothesis  tests  and 
confidence  intervals.  The  validity  of  these  statistics  depends  upon  the  mean  and  variance- 
covariance  model  selected.  The  independent  variable  was  Re,x  with  three  dependent 
variables  (ensemble,  protocol,  and  participant).  For  the  PM  model,  participant  was  set  as 
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the  random-effect  parameter,  and  for  three-way  ANOVA  participant  was  part  of  the  class 
statement. 

Once  significance  was  detected,  Tukey’s  HSD  test  was  used  to  test  all  pair-wise 
comparisons  among  Re  t  means  to  determine  which  ensembles  were  significantly 
different.  Interaction  between  two  variables  was  also  evaluated  (ensemble  x  protocol). 
Significance  levels  were  set  at  a  =  0.05.  Three  hypothesis’  were  tested:  (1)  there  are  no 
differences  between  mean  R0,t  values  among  ensembles,  (2)  there  are  no  differences 
between  mean  Rej  values  among  environments  and  metabolic  rates/demands,  and  (3) 
there  are  no  differences  between  mean  R^j  values  while  AP  changes  within  environments 
and  metabolic  demands. 
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RESULTS 


The  primary  purpose  of  this  paper  was  to  explore  the  methodology  for 
distinguishing  between  garments  based  on  the  total  evaporative  resistance  properties 
within  different  environments  and  work  demands.  The  secondary  purpose  is  to  challenge 
the  relationship  of  Rej  with  respect  to  changes  in  AP  (PSk  -  Pa)  and  Emax.  Experimental 
trials  were  conducted  to  determine  the  evaporative  resistance  for  five  clothing  ensembles. 
The  protocols  included  a  fixed  metabolic  demand  under  three  different  relative  humidity 
levels  for  Phase  1,  and  for  three  metabolic  demands  with  a  fixed  relative  humidity  level 
for  Phase  2.  The  hypothesis’  tested  include:  (1)  there  are  no  differences  between  mean 
Re)T  values  among  ensembles,  (2)  there  are  no  differences  between  mean  R^t  values 
among  environments  and  metabolic  rates/demands,  and  (3)  there  are  no  differences 
between  mean  ReT  values  while  AP  changes  within  environments  and  metabolic 
demands. 

Experimental  Data 

At  the  critical  conditions  measured  data  captured  included  heart  rate  (HR),  rectal 
temperature  (Tre),  skin  temperatures  (calf,  thigh,  upper  arm,  and  chest),  and 
environmental  conditions  (humidity  and  dry,  wet,  and  black  bulb  temperatures).  Oxygen 
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(O2)  consumption  was  measured  and  recorded  at  30  minute  intervals  (at  30,  60  and  90 
minute  point). 

Using  the  measured  data,  other  key  components  were  computed.  The  metabolic 
rate  (M)  was  computed  using  Equation  14  and  the  differences  in  temperatures  and  partial 
pressures  (AT  and  AP)  were  calculated  using  Equations  27-30.  The  data  for  the  critical 
conditions  for  all  the  trials  are  provided  in  Appendix  B  for  Phase  1  and  Appendix  C  for 
Phase  2. 


Total  Insulation 

Results  were  reported  as  Ij.stat  values  and  then  converted  to  Ij ,dyn  by  adjusting  for 
wind  and  movement  as  suggested  by  Havenith  et  al  [32]  as  shown  in  Equation  20. 
Results  for  both  static  and  dynamic  values  are  presented  in  Table  3. 


Table  3.  Total  Insulation  Values  for  Ensembles. 


Clothing  Item 

■^T,stat 

80  W/m2 

Jr.dyn 
160  W/m2 

240  W/m2 

Ensemble  A 

~(WC)~ 

Work  Clothes 

0.180 

0.168 

wmm 

Ensemble  B 

-(CC)- 

Cotton  Coverall 

0.196 

0.182 

0.160 

warn 

Ensemble  C 

~  (PB)  — 

Tyvek  1424 

0.191 

0.178 

0.156 

0.141 

-(PB)- 

Tyvek  1427 

0.190 

0.177 

0.155 

0.140 

Ensemble  D 

-  (WB)  - 

Water  Barrier 

0.189 

0.176 

0.154 

0.140 

Ensemble  E 

~(VB)- 

Vapor  Barrier 

0.185 

0.172 

0.151 

0.137 
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Phase  1 


In  Phase  1,  the  primary  focus  was  to  determine  if  the  methodologies  used  can 
distinguish  differences  among  the  five  selected  ensembles  (WC  -  work  clothes,  CC  - 
cotton  coveralls,  PB  -  particle  barrier,  WB  -  water  barrier,  and  VB  -  vapor  barrier)  and 
evaluate  how  the  environment  affects  Re,T.  There  were  three  different  environments  (R2 
-  hot/dry,  R5  -  moderate,  and  R7  -  warm/humid)  with  a  fixed  moderate  metabolic 
demand  (M2)  of  160  W/m2.  The  average  Rej  values  and  standard  deviations  are 
presented  in  Table  4. 


Table  4.  Phase  1  -  Mean  Re.T  Values  with  Standard  Deviations. 


Ensemble 

All  Data 

Mean  Std  Dev 

R2 

Mean  Std  Dev 

R5 

Mean  Std  Dev 

R7 

Mean  Std  Dev 

A 

0.013 

0.0040 

0.017 

0.0035 

0.012 

0.0029 

0.011 

0.0029 

B 

0.014 

0.0047 

0.018 

0.0046 

0.012 

0.0035 

0.012 

0.0035 

C 

0.015 

0.0052 

0.020 

0.0042 

0.014 

0.0043 

0.013 

0.0047 

D 

0.017 

0.0053 

0.021 

0.0039 

0.016 

0.0051 

0.014 

0.0046 

E 

0.027 

0.0089 

0.034 

0.0100 

0.026 

0.0051 

0.021 

0.0065 

A  Sharpio-Wilkes  statistical  test  for  fit  was  performed  on  all  datasets  to  determine 

the  best  fit  of  the  data  (normal  or  log  normal).  All  of  the  data  fit  well  as  being  normally 

distributed.  The  data  were  analyzed  using  the  mixed  procedure  and  using  a  three-way 

analysis  of  variance  (ANOVA).  The  main  effects  included  three  protocols,  five 

ensembles  and  14  participants.  Not  all  participants  completed  all  trials  and  some  trials 

were  repeated  which  resulted  in  an  unbalanced  design.  Using  SAS  8.1,  the  Mixed  and 

GLM  models  were  used  to  determine  statistical  differences  for  R«  t  among  ensembles, 
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environments,  and  participants.  Participants  were  treated  as  a  blocking  variable.  The 
SAS  code  used  and  data  output  for  Phase  1  is  provided  in  Appendix  D. 

Very  significant  differences  (p<0.0001)  were  found  for  ensemble,  environment, 
and  participant.  Tukey’s  HSD  test  was  performed  to  determine  which  pairs  were 
significantly  different  among  ensembles  and  environments.  This  resulted  with  Ensemble 
E  being  different  from  all  other  Ensembles  and  Ensemble  D  being  different  from 
Ensembles  A  and  B.  This  is  depicted  graphically  along  with  the  mean  R^t  values  for 
ensembles  in  Figure  2.  The  lines  below  the  ensembles  indicate  ensembles  that  are 
statistically  similar. 

Ensemble  E  was  very  different  from  the  other  ensembles  and  could  have 
interfered  with  the  ability  to  differentiate  differences  among  the  ensembles.  Therefore 
the  data  were  analyzed  again  with  Ensemble  E  excluded.  In  this  analysis,  Ensemble  D 
was  different  from  all  other  ensembles  and  Ensemble  A  was  statistically  different  from 
Ensemble  C.  These  data  are  presented  in  Figure  3. 
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Ensemble 


Figure  2.  Phase  1  -  Mean  Re/rby  Ensemble. 


Figure  3.  Phase  1  -  Mean  Re,i  by  Ensemble  w/o  Ensemble  E. 


44 


Interaction  between  ensemble  and  environment  was  tested  and  found  to  be 
significant  (p=0.0187)  with  Ensemble  E  in  the  mix  and  not  significant  (p=0.8820)  when 
analyzed  without  Ensemble  E.  The  interaction  between  environment  and  ensemble  can 
be  seen  graphically  in  Figure  4. 


Figure  4.  Phase  1 — Mean  Re  T  Values  by  Ensemble  and  Environment. 


The  statistical  software  JMP-IN  5.1  was  used  to  analyze  the  mean  Rej  values  for 
each  environment  within  an  ensemble.  This  resulted  with  all  ensembles  being  very 
significantly  different  within  each  environment  (p<  0.001). 


45 


Analysis  of  Rej  by  environment  resulted  with  the  environment  being  significantly 
different.  Tukey’s  HSD  detected  all  pairs  to  be  very  significantly  different  (p<0.001). 
Figure  5  presents  the  mean  Re>T  values  by  environment  and  indicates  that  all 
environments  are  statistically  different  from  the  others. 


R2  R5  R7 


Environment 


Figure  5.  Phase  1 — Mean  Revalues  by  Environment. 


The  statistical  software  JMP-IN  5.1  was  used  to  analyze  the  mean  Re(T  values  for 
each  environment  within  an  ensemble.  This  resulted  the  Rej  values  being  significantly 
different  between  environments  for  all  Ensembles  (p<0.001). 
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Phase  2 


In  Phase  2,  the  primary  focus  was  to  verify  the  methodologies  used  can 
distinguish  differences  among  the  five  selected  ensembles  (WC  -  work  clothes,  CC  - 
cotton  coveralls,  PB  —  particle  barrier,  WB  -  water  barrier,  and  VB  -  vapor  barrier)  and 
evaluate  how  the  metabolic  rate  affects  Re  t-  There  were  three  different  metabolic  rates 
(Ml  -  light  work,  M2  -  moderate  work,  and  M3  -  heavy  work)  with  a  fixed  moderate 
environment  (R5)  at  50%  rh.  The  average  Re, t  values  and  standard  deviations  are 
presented  in  Table  5. 


Table  5.  Phase  2  -  Mean  Re.x  Values  with  Standard  Deviations. 


Ensemble 

All  Data 

Mean  Std  Dev 

Ml 

Mean  Std  Dev 

M2 

Mean  Std  Dev 

M3 

Mean  Std  Dev 

A 

0.011 

0.002 

0.011 

0.002 

0.013 

0.003 

0.011 

0.001 

B 

0.012 

0.003 

0.014 

0.003 

0.012 

0.002 

0.011 

0.003 

C 

0.013 

0.003 

0.015 

0.004 

0.012 

0.002 

0.011 

0.001 

D 

0.015 

0.004 

0.018 

0.005 

0.015 

0.002 

0.012 

0.002 

E 

0.024 

0.006 

0.028 

0.005 

0.024 

0.004 

0.019 

0.003 

The  data  were  reviewed  for  outliers  and  20  out  of  226  data  points  exceeded  the 
mean  ±  two  times  the  standard  deviation.  A  Sharpio- Wilkes  statistical  test  for  fit  was 
performed  on  all  datasets  to  determine  the  best  fit  of  the  data  (normal  or  log  normal).  All 
of  the  data  fit  well  as  being  normally  distributed.  The  data  were  analyzed  using  the 
mixed  procedure  and  using  a  three-way  analysis  of  variance  (ANOVA).  The  main 
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effects  included  three  protocols,  five  ensembles  and  15  participants.  Not  all  participants 
completed  all  trials  and  some  trials  were  repeated  which  resulted  in  an  unbalanced 
design.  Using  SAS  8.1,  the  Mixed  and  GLM  models  were  used  to  determine  statistical 
differences  for  Re/r  among  ensembles,  metabolic  rates,  and  participants.  Participants 
were  treated  as  a  blocking  variable.  The  SAS  code  used  and  data  output  for  Phase  2  is 
provided  in  Appendix  E, 

Very  significant  differences  (p<0.0001)  were  found  for  ensemble,  metabolic  rate, 
and  participant  (p<0.0001).  Tukey’s  HSD  test  was  performed  to  determine  which  pairs 
were  significantly  different  among  ensembles  and  environments.  This  resulted  with 
Ensembles  D  and  E  being  different  from  all  other  ensembles.  There  was  no  statistical 
difference  detected  when  analyzing  the  data  with  and  without  outliers,  therefore  the 
complete  dataset  was  used  for  all  data  references.  Figure  6  depicts  the  mean  Re,x  values 
by  ensembles  graphically.  The  lines  below  the  ensembles  indicate  ensembles  that  are 
statistically  similar. 

Analysis  of  Re/r  by  metabolic  rate  resulted  with  the  environment  being 
significantly  different.  Tukey’s  HSD  detected  all  pairs  to  be  very  significantly  different 
(p<0.0002).  Figure  7  presents  the  mean  Rej  values  by  environment  and  indicates  that  all 
environments  are  statistically  different  from  the  others. 
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Figure  6.  Phase  2  -  Mean  Re/r  by  Ensemble. 


Figure  7.  Phase  2  -  Mean  Re,T  by  Metabolic  Rate. 
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All  trials  were  designed  for  the  participants  to  elicit  a  desired  metabolic  rate  based 
on  varying  the  speed  and  grade  of  the  treadmill  for  each  individual.  During  Phase  2  there 
were  three  desired  metabolic  rates  80  W/m2,  160  W/m2,  and  240  W/m2,  The  average 
metabolic  rates  by  protocol  are  provided  in  Table  6  and  shown  graphically  in  Figure  8. 


Table  6.  Phase  2  -  Average  Metabolic  Rates. 


Ensemble 

Ml 

M2 

M3 

Avg 

A 

121 

175 

250 

183 

B 

118 

177 

241 

OO 

C 

108 

178 

251 

177 

D' 

111 

177 

259 

182 

E 

114 

176 

249 

181 

Average 

114 

176 

250 

Interaction  between  ensemble  and  metabolic  rate  was  tested  and  found  to  be  very 
significant  (p<0,0001).  The  interaction  between  metabolic  rate  and  ensemble  can  be  seen 
graphically  in  Figure  9. 
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Figure  8.  Phase  2  -  Average  Metabolic  Rates  by  Ensemble. 


Figure  9.  Phase  2  -  Mean  Rejby  Metabolic  Rate  and  Ensemble. 
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Looking  for  where  the  interaction  might  occur,  the  statistical  software  JMP-IN 
5.1  was  used  to  analyze  the  mean  Rej  values  for  each  metabolic  rate  within  an  ensemble. 
This  resulted  the  R^t  values  not  being  significantly  different  between  metabolic  rates 
within  Ensemble  A  (p=0.0717)  and  B  (p=0.0610),  and  very  significantly  different  for 
Ensembles  C,  D,  and  E  (p<0.001).  The  complete  JMP-IN  analysis  of  Phase  1  and  2 
(protocols)  is  provided  in  Appendix  F. 
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DISCUSSION 


The  focus  of  this  study  was  to  conduct  experimental  trials  to  explore  two  research 
areas.  First,  trials  were  conducted  to  distinguish  among  garments  based  on  the  total 
evaporative  resistance  properties  between  different  environments  and  work  demands. 
Experimental  trials  were  conducted  in  each  phase  to  determine  the  evaporative  resistance 
for  five  selected  clothing  ensembles.  The  protocols  included  a  fixed  metabolic  demand 
under  three  different  relative  humidity  levels  for  Phase  1,  and  for  three  metabolic 
demands  with  a  fixed  relative  humidity  level  for  Phase  2.  Second,  the  data  from  the 
experimental  trials  were  used  to  discern  weather  or  not  the  generally  accepted  theory  that 
Re,T  remains  constant. 

Internal  Validity 

Phase  1  and  2  both  had  one  protocol  that  had  the  same  design  —  a  moderate  work 
rate  (M2  -  160  W/m2)  and  a  moderate  environment  (R5  -  50%  rh).  Ensemble  C  (PB) 
was  changed  between  Phase  1  and  2,  but  the  other  ensembles  remained  the  same.  Some 
of  the  same  participants  from  Phase  1  were  used  again  in  Phase  2,  but  most  were 
different.  Comparing  the  moderate  work  rate  and  moderate  environment  (M2R5)  data 
from  both  phases  provided  internal  validity  to  the  methodology.  The  mean  values  for 
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ensembles  for  Phase  1  and  2  were  plotted  and  presented  in  Figure  10,  The  average 
metabolic  rates  for  Phase  1  and  2  are  provided  in  Figure  11. 


0.030 

0.025 

0.020 
£  0.015 
0.010 

0.005 

0.000 

A  B  C  D  E 

Ensemble 

Figure  10.  Comparison  of  M2R5  Mean  Rej. 


In  Figure  10  the  grouping  of  data  points  for  each  ensemble  appears  to  be  tightly 
correlated  with  the  exception  of  Ensemble  C.  As  discussed  previously,  Ensemble  C  was 
changed  from  a  Tyvek  1424  for  Phase  1  to  Tyvek  1427  for  Phase  2.  Statistical  analysis 
using  JMP-IN  5.1  was  used  to  analyze  the  M2R5  data  to  compare  mean  Re>x  values 
within  ensembles.  This  analysis  resulted  with  only  Ensemble  C  being  significantly 
different  (p=0.0349).  The  statistical  results  (p  values)  and  mean  Rej  values  are  provided 
in  Table  8.  The  JMP-IN  analysis  is  provided  in  Appendix  G. 
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Figure  1 1 .  Comparison  of  M2R5  Mean  Metabolic  Demands. 


Table  7.  Statistical  Analysis  of  M2R5  Dataset. 


Ensemble 

p  Value 

Mean  ReX 

Phase!  Phase  2 

A 

0.7692 

B 

0.4741 

0.013 

0.012 

C 

0.0349 

0.015 

0.012 

D 

0.2577 

0.017 

0.015 

E 

0.1653 

0.026 

0.024 

At  first  glance,  the  metabolic  rates  seen  in  Figure  1 1  appear  to  be  significantly 
different.  However,  the  scale  is  compressed  making  the  small  differences  (<  10%)  seem 
larger.  The  differences  in  Phase  1  and  2  were  not  enough  to  change  the  conformation  of 
internal  validity. 
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Comparison  to  Other  Studies 


Total  Insulation 

In  order  to  determine  the  evaporative  resistance,  an  understanding  of  the  ensemble 
properties  must  be  understood.  The  clothing  properties  were  derived  from  manikin 
experiments  conducted  at  the  Institute  for  Environmental  Research,  Kansas  State 
University  by  Dr.  Elizabeth  McCullough.  Using  her  manikin,  and  following  ASTM  F 
1291,  she  was  able  to  determine  the  total  clothing  insulation  (Ix.stat),  intrinsic  clothing 
insulation  (Ici)  and  the  clothing  area  factor  (fci)  for  the  six  ensembles  used  in  the 
experimental  trials.  As  reported  by  Havenith  et  al.  [25]  heated  manikin  results  for 
standing/no  wind  appears  to  be  on  average  0.023  °C  m2/W  higher  than  human 
participants.  While  some  studies  support  this  claim,  other  studies  find  the  manikin  data 
as  being  lower  than  human  participants.  After  adjusting  for  wind  and  movement  the  Ij.dyn 
values  were  compared  to  other  studies  that  had  similar  ensembles.  The  studies  included 
Barker  et  al.  [23],  Kenney  et  al.  [21],  and  Bernard  and  Matheen  [9].  After  adjusting  the 
It  values,  the  values  used  in  this  study  are  clearly  higher  than  other  studies  for  similar 
ensembles.  The  reported  It  values  are  shown  in  Table  9. 
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Table  8.  It  Values  from  Different  Studies. 


IT 

(m2  K/W) 

Ensemble 

Current 

Barker 

Kenney 

Bernard 

A 

0.147 

0.084 

0.050 

B 

0.160 

0.107 

0.056 

0.107 

C 

0.156 

0.086 

0.059 

D 

0.154 

0.086 

0.050 

E 

0.151 

0.086 

0.035 

The  Barker  et  al.  and  Bernard  and  Matheen  studies  reported  It  values  that  were  33 
-  45%  lower,  and  the  Kenney  et  al.  study  reported  values  that  were  62  -  77%  lower.  A 
primary  difference  in  all  of  these  studies  is  the  adjustment  for  wettedness.  Although 
there  isn’t  a  set  standard  for  adjusting  for  clothing  wettedness  thus  far,  many  researchers 
use  a  50%  default  adjustment.  The  Barker  at  al.  and  Bernard  and  Matheen  studies  both 
used  a  45%  adjustment  for  wettedness.  Had  the  current  values  been  adjusted  for 
wettedness,  the  It  values  would  match  the  Barker  et  al.  study  very  well.  On  the  other 
hand,  Kenney  et  al.  used  a  simultaneous  derivation  method  to  compute  IT.  Using  his 
methodology  to  compute  It  with  the  Phase  1  data  resulted  in  too  much  variation  in  It  to 
make  it  useful. 


On  the  face  of  it,  having  a  good  estimate  of  the  It  is  important  because  it  is  used 

to  compute  Re/r-  However,  Barker  et  al.  demonstrated  that  relatively  large  changes  in  It 

result  in  minor  changes  in  R^t-  Therefore  the  manikin  data  adjusted  for  wind  and 

movement  is  sufficient  for  determining  the  It  of  the  ensembles  used. 
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Total  Evaporative  Resistance 

In  1993  Kenney  at  al.  [19],  building  from  previous  research,  setup  the  framework 
for  conducting  human  experiments  in  a  climate  controlled  heated  chamber.  The 
methodology  they  established  is  still  in  use. 

Using  the  principle  of  the  prescriptive  zone  as  established  by  Lind  [15],  the 
determination  of  the  inflection  point  is  established  by  selecting  the  point  preceding  a  rise 
in  Tre.  At  the  inflection  point,  critical  conditions  exist  where  S  =  0  and  E^x  =  Ereq. 

From  these  conditions  the  basic  heat  balance  equation  can  be  manipulated  by  substituting 
terms  for  Emax  and  Ereq  and  solving  for  Rej. 

Similarly  to  the  It  values,  mean  R^j  values  were  compared  to  the  same  studies  - 
Barker  et  al.  [23],  Kenney  et  al.  [21],  and  Bernard  and  Matheen  [9],  Although  there 
were  large  differences  (33  -  45%  lower  It  values  reported  by  Barker  et  al.  and  62  -  77% 
lower  values  reported  by  Kenney  et  al.),  the  Rej  values  had  less  differences  as  shown  in 
Table  10. 
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Table  9.  Re.T  Values  from  Different  Studies. 


Re,T 

(kPa  m2/W) 

Ensembl 

e  Current  ; 

Barker 

Kenney 

Bernard 

A 

0.0133 

0.0131 

0.0092 

B 

0.0140 

0.0159 

0.0096 

0.0155 

C 

0.0154 

0.0163 

0.0112 

D 

0.0174 

0.0176 

0.0123 

E 

0.0273 

0.0136 

0.0344 

Barker  et  al.  reported  three  Re, t  values  that  were  within  6%,  one  at  14%  and  one  at 
50%  (vapor  barrier  suit).  Kenney  et  al.’s  values  ranged  from  26  -  32%  difference,  and 
Bernard  and  Matheen’s  reported  value  was  1 1%  higher.  Barker’s  Ij  values  were  close  to 
the  ensembles  used  in  this  study  with  the  exception  of  not  adjusting  these  values  by  45% 
for  account  for  wettedness.  However,  even  with  the  45%  difference  in  IT,  there  is  minor 
differences  in  R^t-  Barker  et  al.  had  previously  reported  this  relationship,  and  this  study 
supports  it.  Based  on  the  Barker  et  al.  and  Bernard  and  Matheen  studies,  the  Re, T  values 
calculated  in  this  study  appear  in  line  with  other  research. 

Phase  1 

The  methodology  used  was  able  to  distinguish  effectively  between  the  selected 

ensembles  as  illustrated  in  Figures  2  and  3  showing  significant  differences  among 

Ensembles  D  and  E.  Since  Ensemble  E  is  a  vapor  barrier  suit,  it  is  expected  to  be 

different  from  ensembles  that  do  not  prevent  vapor  transmission  such  as  cotton  and 

Tyvek  coveralls.  Similarly,  Ensemble  D  is  a  liquid  barrier  suit,  so  it  is  also  expected  to 

59 


be  different  from  particle  barrier  and  cotton  clothing  (Ensembles  A,  B  and  C )  with 
respect  to  evaporative  resistance. 

In  Figure  4  there  is  a  clear  difference  between  environments.  The  differences 
between  the  mean  Re  x  values  remains  the  same  within  Ensembles  A  -  D,  but  increases 
for  Ensemble  E.  The  increase  difference  accounts  for  the  interaction  between  the 
environment  and  ensemble  and  is  verified  by  not  seeing  an  interaction  when  the  data  is 
analyzed  without  Ensemble  E. 

The  fact  that  there  is  a  difference  between  environments  is  by  itself  an  important 
finding.  The  relationship  between  Re,T  and  AP  with  respect  to  Emax  (Equation  1 0)  has 
generally  been  accepted  that  Re,T  remains  constant  as  AP  and  Emax  change.  This 
relationship  is  alluded  to  in  ISO  7933  and  discussed  by  Parsons  (2003)  [7]. 

Again,  Figure  4  clearly  shows  that  RCjt  is  not  the  same  as  the  environment 
changes.  Using  the  data  from  Phase  1  and  Phase  2,  AP  was  plotted  against  Rejto  test  this 
theory.  The  data  were  plotted  for  each  ensemble  across  all  protocols  and  a  regression  line 
was  calculated.  All  of  the  graphs  are  presented  in  Appendix  H,  and  the  graph  for  Phase  1 
Ensemble  A  is  shown  in  Figure  12.  It  is  obvious  that  Re,T  does  not  remain  constant  as  AP 
changes,  and  the  rate  of  change  (slope)  appears  constant  within  Phase  1  for  the  different 
ensembles  suggesting  that  environment  is  a  factor.  However,  the  Phase  2  graphs  do  not 
show  any  consistent  effect  for  activity.  The  Phase  2  data  is  confounded  by  the  metabolic 
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rate  and  therefore  there  isn’t  an  expected  effect.  Regression  analysis  performed  on  the 
data  resulted  with  Phase  2  ensembles  having  an  R-square  value  of  0.002  -  0.1 15,  while 
the  Phase  1  ensembles  ranged  from  0.388  -  0.617. 


Figure  12.  Effect  of  Environment:  Re  ?  vs.  AP  -  Ensemble  A. 


Table  10  presents  the  slopes  and  intercepts  of  the  regression  analysis.  By 
reviewing  the  slope  values,  environment  (Phase  1)  appears  to  have  consistent  and 
significant  slope  values  across  the  ensembles.  However,  the  values  for  Phase  2  suggest 
the  slopes  are  neither  consistent  nor  significant.  Although  Ensemble  A  in  Phase  2 
presents  a  good  positive  slope,  all  of  the  other  ensembles  do  not. 
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Table  10,  Regression  Analysis  -  AP  by  Re,x 


Ensemble 

R2 

Phase  1 
Slope 

Intercept 

R2 

Phase  2 
Slope 

Intercept 

A 

0.617 

0.0001 

0.0004 

0.115 

0.0244 

0.0001 

B 

0.453 

0.0001 

0.0025 

0.018 

0.3951 

0.0007 

C 

0.477 

0.0001 

0.0037 

0.001 

0.8478 

0.0001 

D 

0.388 

0.0001 

0.0002 

0.001 

0.8420 

0.0009 

E 

0.534 

0.0001 

0.7474 

0.005 

0.6430 

0.0006 

Phase  2 

Similar  to  Phase  1,  the  methodology  used  was  able  to  distinguish  effectively 
between  the  selected  ensembles  as  illustrated  in  Figure  6,  showing  very  significant 
differences  among  Ensembles  D  and  E  as  compared  to  Ensembles  A,  B,  and  C.  Figure  7 
indicates  Rg/r  decreases  as  the  metabolic  rate  increases.  The  interaction  between 
ensemble  and  metabolic  rate  is  clearly  seen  by  observing  the  differences  between 
metabolic  rates  within  each  ensemble  increase  corresponding  to  the  reduction  to  Rej. 

There  is  not  much  difference  in  the  mean  Revalues  between  the  metabolic  rates 
within  Ensemble  A  indicating  good  evaporative  cooling  (high  permeability).  However, 
progressing  through  the  ensembles,  the  differences  between  protocols  increases 
indicating  the  metabolic  rate  (activity)  plays  an  increasing  role  in  lowering  the  mean  Re/r. 
Again,  this  relationship  was  verified  by  using  JMP-IN  to  test  the  differences  between  the 
protocols  for  each  ensemble.  Ensemble  A  and  B  were  not  significantly  different  whereas 
the  others  (Ensembles  C,  D,  and  E)  all  resulted  as  being  very  significantly  different. 
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As  discussed  previously,  increasing  the  metabolic  rate  results  with  decreasing  the 
mean  Rej  values.  This  effect  is  presented  differently  in  Figure  13  where  the  difference 
among  ensembles  is  distinct  for  Ml,  but  changes  as  the  metabolic  demand  increases.  For 
Ml,  there  appears  to  be  a  step  effect  between  the  ensembles.  In  M2  and  more  so  in  M3, 
this  step  effect  disappears  as  the  ensembles  appear  to  reach  the  lower  limit  of  Rej. 


Figure  13.  Effect  of  Metabolic  Rate  on  Re/r. 


Havenith  et  al.  [30]  and  Parsons  et  al.  [31]  explain  the  relative  decrease  in  Re  t  due 
to  the  increased  air  movement  through  the  clothing.  They  use  the  term  pumping  action  to 
explain  that  as  an  individual  moves,  air  is  pumped  into  and  out  of  their  clothing.  As  the 
air  moves  through  the  clothing,  the  effective  It  and  Rej  decreases  due  to  increased 
convective  and  evaporative  cooling. 
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CONCLUSIONS 


The  primary  purpose  of  this  research  was  to  explore  the  differences  among 
garments  based  on  the  total  evaporative  resistance  properties  among  different 
environments  and  work  demands.  The  secondary  purpose  was  to  challenge  the 
relationship  of  Rej  with  respect  to  changes  in  AP  and  Emax.  Experimental  trials  were 
conducted  to  determine  the  evaporative  resistance  for  five  clothing  ensembles  per  phase. 
The  protocols  included  a  fixed  metabolic  demand  under  three  different  environments 
(levels  of  relative  humidity)  for  Phase  1,  and  a  fixed  relative  humidity  level  with  three 
metabolic  demands  with  for  Phase  2.  The  fundamental  step  in  these  studies  was  being 
able  to  distinguish  the  point  just  before  the  transition  of  compensable  heat  stress  to 
uncompensable  heat  stress  (Ereq  =  Emax). 

Statistical  analysis  of  the  data  showed  that  the  methodology  used  was  able  to 
distinguish  well  among  the  selected  ensembles.  Data  from  Phase  1  found  that  Ensemble 
E  was  different  from  all  others  and  Ensemble  D  was  different  from  A  and  B,  More 
importantly,  the  data  revealed  a  relationship  with  the  environment.  The  mean  Revalues 
for  each  ensemble  decreases  as  the  humidity  increased.  The  changes  to  Rc  j  due  to 
environment  were  explored  further. 
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The  default  assumption  has  been  that  Rej  remains  constant  as  AP  changes.  This 
relationship  between  Re/r  and  AP  was  challenged  and  found  that  Re>T  does  not  stay 
constant  as  generally  accepted.  Environment  (relative  humidity)  effects  Re/r  as  well  as 
AP.  This  relationship  needs  to  be  studied  further  before  it  is  fully  understood. 

The  Phase  2  analysis  resulted  with  Ensembles  D  and  E  being  different  from  all 
other  ensembles.  As  expected,  with  increased  activity  mean  R^t  values  decreased. 
Ensembles  D  and  E  had  the  biggest  decreases  in  Re/r,  while  Ensembles  A,  B,  and  C 
appeared  to  reach  a  lower  limit  associated  with  the  ensemble  permeability  properties. 
The  decrease  in  Re/r  from  metabolic  demand  was  related  to  the  pumping  action  of  air 
through  the  ensemble  from  movement. 

The  null  hypothesis’  were  rejected  for  all  three  hypothesis’  tested.  The  data 
shows  (1)  there  are  differences  between  Re  j  values  within  ensembles,  (2)  there  are 
differences  between  Re/r  values  within  ensembles  and  between  the  different  metabolic 
rates/demands,  and  (3)  Re/r  does  not  remain  constant  while  AP  changes. 
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Table  Al.  Characteristics  of  Participants  in  Experimental  Trials. 


Height  (cm) 

Weight  (kg) 

Surface 
Area  (m2) 

Y1S1 

M 

26 

180 

95 

2.14 

Y1S2 

F 

26 

163 

52 

1.55 

Y1S3 

M 

24 

183 

86 

2.08 

Y1S4 

M 

25 

183 

77 

1.99 

Y1S5 

F 

23 

152 

63 

1.59 

Y1S6 

F 

27 

170 

91 

2.02 

Y1S7 

M 

35 

189 

101 

2.28 

Y1S8 

F 

39 

155 

46 

1.42 

Y1S9 

M 

20 

183 

130 

2.48 

Y1S10 

M 

30 

191 

110 

0.00 

Y1S11 

M 

32 

173 

71 

1.84 

Y1S12 

M 

43 

178 

112 

2.28 

Y1S13 

M 

28 

185 

95 

2.19 

Y1S14 

F 

44 

165 

65 

1.72 

Y2S1 

F 

27 

163 

52 

1.55 

Y2S2 

M 

28 

185 

95 

2.19 

Y2S3 

F 

27 

170 

91 

2.02 

Y2S4 

M 

26 

180 

95 

2.15 

Y2S5 

M 

27 

175 

98 

2.13 

Y2S6 

M 

20 

180 

83 

2.03 

Y2S7 

M 

20 

183 

72 

1.93 

Y2S8 

M 

24 

163 

64 

1.68 

Y2S9 

M 

43 

149 

75 

1.69 

Y2S10 

M 

49 

175 

86 

2.02 

Y2S11 

F 

18 

170 

56.8 

1.66 

Y2S12 

F 

20 

157 

56.8 

1.57 

Y2S13 

M 

21 

185 

81.8 

2.06 

Y2S14 

M 

22 

175 

66 

1.80 

Y2S15 

M 

28 

185 

86 

2.11 

Average  ±  Std  Dev 

28.3  ±8.1 

173.9  ±11.5 

81.1  ±20.1 

1.87  ±0.45 
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Appendix  B 

Experimental  Data  -  Phase  1:  Data  Dictionary 

Title  Description 

Code  Participant  Code 

Gender  Gender  of  participant 

Protocol  Design: 

Proto  Environment  (R2  (20%  rh),  R5  (50%  rh),  R7  (70%  rh))  or 

Metabolic  Demand  (Ml  (80  W/m2),  M2  (160  W/m2),  M3  (240  W/m2)) 
gns  Ensemble:  (A  (work  clothes),  B  (cotton  coveralls),  C  (particle  barrier), 

D  (liquid  barrier),  E  (vapor  barrier)) 

Tdb  Ambient  air  temperature  (dry  bulb)  in  degrees  Celsius 

Tpwb  Wet  bulb  air  temperature  in  degrees  Celsius 

Tg  Black  bulb  air  temperature  in  degrees  Celsius 

S(m/s)  Speed  in  meters  per  second 

G(%)  Grade  of  treadmill  in  percentage 

HR  Heart  rate 

Tre  Body  core  temperature  (rectal) 

Tch  Skin  temperature  at  the  chest 

Tarm  Skin  temperature  at  the  upper  arm 

Tth  Skin  temperature  at  the  thigh 

Tcalf  Skin  temperature  at  the  calf 

Met  Calculated  metabolic  work  based  on  02  consumption  in  Watts 

BSA  Body  surface  area  in  square  meters 

MSA  Met  divided  by  the  BSA  (W/m2) 

Tsk  Average  Skin  temperature 

Psk  Partial  pressure  of  the  water  vapor  at  the  skin 

Pv  Partial  pressure  of  the  water  vapor  in  the  air 

Psk-Pv  AP:  Difference  between  Psk  and  Pv 
Tair-Tsk  AT:  Difference  between  Tdb  and  Tsk 

ReT  Total  evaporative  resistance 
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Experimental  Data  -  Phase  1 
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Experimental  Data  -  Phase  1 
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Experimental  Data  -  Phase  1 
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Experimental  Data  -  Phase  1 


31.9  22.8  30.S  1.03  0  89  37.4  36.6  36.3  36.3  36.2  36.4  6.1  2.2  2.38 


Code  Gender  Proto  Ens  |  Tdb  Tpwb  Tg  S(m/s)  G(%)  HR  Tre  Tch  Tarm  Tth  Tcalf  |  Tsk  Psk  Py  BSA  Met  MSA  Psk-Pair  Tair-Tsk  ReT 


Sll  M  R5  E  30.9  23.5  29.8  1.02  0  117  37.9  36.9  36.9  36.6  36.7  36.8  6.2  2.4  2.28  415  182  3.8  -5.9  0.026 
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Appendix  C 

Experimental  Data  -  Phase  2:  Data  Dictionary 


Title 

Description 

Code 

Participant  Code 

Gender 

Gender  of  participant 

Proto 

Protocol  Design; 

Metabolic  Demand  (Ml  (80  W/m2),  M2  (160  W/m2),  M3  (240  W/m2)) 

Ens 

Ensemble:  (A  (work  clothes),  B  (cotton  coveralls),  C  (particle  barrier), 
D  (liquid  barrier),  E  (vapor  barrier)) 

Tdb 

Ambient  air  temperature  (dry  bulb)  in  degrees  Celsius 

Tpwb 

Wet  bulb  air  temperature  in  degrees  Celsius 

Tg 

Black  bulb  air  temperature  in  degrees  Celsius 

S(m/s) 

Speed  in  meters  per  second 

G(%) 

Grade  of  treadmill  in  percentage 

HR 

Heart  rate 

Tre 

Body  core  temperature  (rectal) 

Tch 

Skin  temperature  at  the  chest 

Tarm 

Skin  temperature  at  the  upper  arm 

Tth 

Skin  temperature  at  the  thigh 

Tcalf 

Skin  temperature  at  the  calf 

Met 

Calculated  metabolic  work  based  on  02  consumption  in  Watts 

BSA 

Body  surface  area  in  square  meters 

MSA 

Met  divided  by  the  BSA  (W/m2) 

Tsk 

Average  Skin  temperature 

Psk 

Partial  pressure  of  the  water  vapor  at  the  skin 

Pv 

Partial  pressure  of  the  water  vapor  in  the  air 

Psk-Pv 

AP:  Difference  between  Psk  and  Pv 

Tair-Tsk 

AT:  Difference  between  Tdb  and  Tsk 

ReT 

Total  evaporative  resistance 
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Code  Gender  Proto  Ens  Tdb  Tpwb  Tg  S(m/s)  G(%)  HR  Trt  Tch  Tarm  Tth  Tcalf  Met  BSA  MSA  Tsk  Psk  Pv  Psk>Pair  Tair-Tsk  R^t 

1  F  Ml  A  44.90  34.00  43.20  0.43  0.0  86  37.06  36.27  35.95  36.56  36.83  170  1.55  110  36.34  6.05  4.59  1.5  8.6  0.009 
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Code  Gender  Proto  Ens  Tdb  Tpwb  Tg  S(m/s)  G(%)  HR  Tre  Tch  Taro  Tth  Tcalf  Met  BSA  MSA  Tsk  Psk  Pv  Fsk-Fair  Tatr-Tsk  R,T 

~3  F  Ml  A  4120  34.50 . 44.00  0M  OO  106  37,57  36.68  3174 . 37.40  ~~36~94"  249  2.02  123  36.89  6.24  4.76  1.5  8.3  0.009 


29.90  22.90  28.50  1.27  0.0  95  37.82  35.76  36.27  32.29  35.51  464  2.16  216  35.17  5.67  2.32  3.4  -5.3 

28.20  20,20  26.70  1.52  1.0  111  37,91  35.43  35.52  35.72  34,74  561  2.16  261  35.38  5.74  1.83  3.9  -7.2 
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36.70  26.00  35.50  1.31  0.0  119  37.39  35.91  35.03  36.07  35,70  270  1.68  161  35.64  5.82  2,64  3.2  1.1  0.019 

35.20  23,80  33,10  1.64  2,0  161  38.44  36.81  37.08  35.77  35,84  408  1.68  243  36.49  6.10  2.18  3.9  4.3  0,017 

33,80  25.00  33.00  0.41  0.0  117  37.64  36.73  37.11  36,74  36.99  206  1.68  123  36.90  6.24  2.58  3.7  -3.1  0,035 

31.00  21.90  28.60  1.32  0.0  140  37.92  35.94  37.88  36.90  36,45  334  1.68  199  36.82  6.21  2.02  4.2  -5.8  0.026 

29.90  21.05  28.05  1.62  2.0  135  38.00  35.39  35.31  36,10  36,57  360  1.68  214  35.74  5.86  1.90  4.0  -5.8  0,023 


Cod?  Gender  Proto  Ens  I  Tdb  Tpwb  Tg  S(m/s)  G  (%)  HR  Tre  Tch  Tarm  Tth  Tcalf  Met  BSA  MSA  Tsk  Psk  Pv  Psk-Pair  Tair-Tsk 


26.60  19.90  27.20  1.55  2.5  140  38.53  36.82  37.45  35.87  36.33  592  2.02  293  36.72  6.18  1.87  4.3  -10.1 


Code  Gender  Proto  Ens  Tdb  Tpwb  Tg  S(m/s)  G(%)  HR  Tre  Tch  Tairo  Tth  Tcalf  I  Met  BSA  MSA  Tsk  Psk  Pv  Psk-Pair  Tair-Tsk 
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SAS  Code  -  Phase  1 

options  nodate  nonumber; 

libname  Vc  1 F: \USF\NTOSH  Studies\evap  res  YrlV; 

*  SAS  Code  for  Analyzing  Re,T  for  Phase  1; 

%macro  meanl  (varl,  var2,  var3,  var4) ; 

Proc  Means  data=Vc,ret  n  mean  var  std  stddev; 
title  M SAS  Analysis  of  Phase  1  Data"; 

Class  £var2  &var3  &var4; 
var  &varl; 

Run; 

%mend; 

%msanl  (ReT,  ensemble)  ; 

%meanl  (ReT,  ensemble,  proto) ; 

%meanl  (ReT,  proto) ; 

%macro  anovl  (varl,  var 2,  var 3,  var4); 

Proc  glm  data=vc.ret; 

title  “Three  way  ANOVA  using  Proc  GLM  for  &varl  Data”; 

Class  &var2  &var3  &var4; 

Model  &varl  =  &var2  &var3  &var4; 

Ismeans  &var2  &var3  &var4  /pdiff  adjust=Tukey  alpha=0*G5; 
run; 

%mend; 

%anovl  (ReT,  ensemble,  proto,  sub j ) ; 

%macro  anov2  (varl,  var2,  var3,  var4); 

Proc  glm  data=vc.ret; 

title  "Three-way  ANOVA  of  &varl  data  set:  Testing  Interaction  of 
&var2  x  &var3”; 

Class  &var2  &var3  &var4; 

Model  &varl  -  &var2  |  &var3  &var4; 

*lsmeans  &var2  I  &var3  /pdiff  adjust-Tukey  alpha™G.05;  - 
run; 

%mend; 

%anov2  (ReT,  ensemble,  proto,  sub j ) ; 

%macro  mixedl  (varl,  var 2,  var 3,  var 4) ; 

Proc  mixed  data=vc.ret; 

title  “Analysis  of  $varl  using  the  Mixed  Model”; 

Class  &var2  &var3  &var4; 

Model  Svarl  =  &var2  &var3; 

Random  &var4; 

LSmeans  &var2  &var3  /adjust=tukey  alpha=,05; 
run; 

%mend; 

%m±xsdl  (ReT,  ensemble,  proto,  sub j ) ; 
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Number  of  observations  220 


Three  way  ANOVA  using  Proc  GLM  for  ReT  Data 
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Three-way  ANOVA  of  ReT  data  set:  Testing  Interaction  of  ensemble  x  proto 
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Ensemble  3  0.00038033  0.00012678  14.49  <.0001 
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Appendix  E 


SAS  Code  -  Phase  2 

options  nodate  nonumber; 

libname  Vc  1 F : \USF\NIOSH  Studies\evap  res  Yr2  V  ; 

*  SAS  Code  for  Analyzing  Re , T  for  Phase  2; 

%macro  meanl  (varl,  var2,  var3,  var4); 

Proc  Means  data=Vc.ret  n  mean  var  std  stddev; 
title  "SAS  Analysis  of  Pase  2  Data"; 

Class  &var2  &var3  &var4; 
var  &varl; 

Run  ; 

%mend; 

%meanl  (ReT,  ensemble) ; 

%meanl  (ReT,  ensemble,  M)  ; 

%  meanl  (ReT,  M)  ; 

%macro  anovl  (varl,  var2,  var3,  var4); 

Proc  glm  data=vc.ret; 

title  "Three  way  ANOVA  using  Proc  GLM  for  &varl  Data"; 

Class  &var2  &var3  &var4; 

Model  &varl  =  &var2  &var3  &var4; 

lsmeans  &var2  &var3  &var4  /pdiff  adjust=Tukey  alpha=0.05; 
run  ; 

%mend; 

%anovl  (ReT,  ensemble,  M,  sub j ) ; 

%macro  anov2  (varl,  var2,  var3,  var4) ; 

Proc  glm  data=vc.ret; 

title  "Three-way  ANOVA  of  &varl  data  set:  Testing  Interaction  of 
&var2  x  &var3” ; 

Class  &var2  &var3  &var4; 

Model  &varl  =  &var2  |  &var3  &var4; 

lsmeans  &var2  |  &var3  /pdiff  adjust=Tukey  alpha=0.05; 
run; 

%mend; 

%anov2  (ReT,  ensemble,  M,  sub j ) ; 

%macro  mixedl  (varl,  var2,  var3,  var4) ; 

Proc  mixed  data=vc.ret; 

title  "Analysis  of  $varl  using  the  Mixed  Model"; 

Class  &var2  &var3  &var4; 

Model  &varl  =  &var2  &var3; 

Random  &var4; 

LSmeans  &var2  &var3  /adjust=tukey  alpha=.05; 
run; 

%mend; 

%mixedl  (ReT,  ensemble,  M,  sub j ) ; 
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Appendix  F 


Re,T  Response  to  Protocol  (Environment)  by  Ensemble 


Root  Moan  Square  Error  0.00326 
Moan  of  Response  0.013052 

Observations  (or  Sum  Wgts}  42 


m 


0.01113333  0.00054176  0.011133 


[  Analysis  of  Variance 


Source  OF  Sum  of  Squares  Mean  Square  F  Ratio 

Model  2  0.00039136  0.00019$  18.4112 

Error  39  0.00041452  0.000011  |Prob»F 

&  Total  41  0.00080590  <.0001 

( Parameter  intimates  ” '"} 

Tenm  Estimate  Std Error  tRefio  Prob>ft| 

Intercept  0,01412$  0.000504  28.03  <0001 

mi*im  0.0042596  0.000728  $.87  <0001 

Pmio|R5|  -0.001288  0,000712  AM  0.0828 

Source  Nparm  DF  Sum  of  Squares  F  Ratio  jprob  >  F 

Proto  2  Z  0.00030138  18.4112  <0001 
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Re/r  Response  to  Protocol  (Environment)  by  Ensemble 


|  Ensemble«C  . 

[Whole  Model 

[  Actual  by  Predicted  Plot 

"I  . . 


Ret  Predicted  P-0.0002  RSq*0.35 

j  Summary  of  Fit  J 

RSquare  •  0.346162 

RSquareAdj  0.315591 

Root  Moan  Square  Error  0.00461 1 

Moan  of  Response  0.015814 

^Observations  (or  Sum  V&!$)  43 

(liw^eis  of  Variance  -  '  .  j 

Source  OF  Sum  of  Squares  Mean  Square  F  Ratio 

Model  2  0,0004  $421  0.000227  10.6834 

Error  40  0,00085030  0 .000021  Pm  >  F 

C,  Total  42  0.00130451  0,0002 

['Parameter  Estimates  TTyT] 


[leverage  Plot 


m 

!  0.026- 

1 0.020- 
$ 

J  0.015  ■ 


^  0.010-1- 


.01 2 .01 3.01 4 .01 6 .016.01 7 .0 1 8 .01 9 .020 
Proto  Leverage,  P*0.ODD2 

[Least  Squares  Means  Table  J 

Level  Least  Sq  Mean  Sid  Error  Mean 
R2  0,02028571  0.00123223  0.020286 


0.01485714 

0,01253333 


0,00123223  0,014857 
0.0011004$  0,012533 


.Term  Estimate  Std  Error  t  Ratio  Prob»jtf 

Intercept  0.0138821  0060703  22.58  <.0001 

Protop2)  0.0043837  O.0O1O01  4.39  <,0001 

ProtofRSJ  *0.00103$  0,001001  *1JS  0.3072 

[Effect  Tests _  ' 

Source  Nparm  OF  Sum  of  Squares  F  Ratio  Prob  >  F 

PmtO  2  2  0.00045421  10.6834  0.0002 
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Appendix  F  (Continued) 


Rej  Response  to  Protocol  (Environment)  by  Ensemble 


RSqoarc  0.315072 

&  Square  Adj  0.28321  $ 

Root  Mem  Square  Enter  0,064662 

Mean  of  Response  0,01 78? 

Observations  (or  Sum  Wgts)  46 


R2  0.0220866?  0,00121020  0.02206? 

R$  0.01688333  0.60110476  0.016833 

R?  0,01446164  0.60126007  0.014462 


|  Analysis  of  Variance  ] 

Source 

OF  Sum  of  Squares 

Mean  Square 

F  Ratio 

Model 

2  0.000434$$ 

0,00021? 

96902 

Error 

43  0.00094466 

0.000022 

Proh  >  F 

O.  Total 

4$  0.00137322 

0.0003 

f  Parameter  Estimates 

) 

Term 

Estimate  Std  Error 

t  Ratio 

PfOfe>ftf 

Intercept 

0.0177872  0.060697 

2S.S1 

<0001 

ProioiR2J 

0.0042795  0.000987 

4,34 

<0001 

Pmtop$I 

*0,000954  0.060945 

*1.01 

0,3184 

[  Effect  Tests  j 

Source  Nparm  OF  Sum  of  Squares 
Rmto  2  2  0.0004348$ 


¥  Ratio  Prob  >  f 
0.8902  0,0003 
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Appendix  F  (Continued) 


Re  j  Response  to  Protocol  (Environment)  by  Ensemble 


[Ensemble*! 

[Response  ReT  . : .  . 1 .  j 

(Whole Model  . ■ 

Proto 

. : . ) 

Actual  by  Predicted  Plot 

{ Leverage  Plot 

. . :  J 

0.06- 


Ret  Predicted  P«.0O01  RSq*0.38 
RMSE*0.O07 


0.06* 


0,O5*H 


.030 

Ptoxo  Leverage.  P<.0O0i 

Least  Squares  Means  Table 


.035 


i.  Summary  of  Fit 

.  \ 

] 

Level 

Least  Sq  Mean 

Sid  Error 

Mean 

RSquare 

0.370396 

R2 

0.03281250 

0.0017687! 

0.032612 

RSquare  Ad) 

0.349843 

R3 

0,02613333 

0.00181038 

0.026133 

Root  Mean  Square  Error 

0.00703$ 

R7 

0.01976571 

0,00188014 

0.019788 

Mean  of  Response 

0.O26533 

Observations  (or  Sum  Vt&ts} 

4$ 

— ! - : — ^ 

1  Analysis  of  Variance 


Source  DF  Sum  of  Squares  Moan  Square  F  Ratio 
Mode!  2  0.00127067  O. 00063$  12,6380 

Error  42  0.00207853  0.0000  40  Pmfc»F 

C,  Total  44  0,00334920  <,0001 


|  Parameter  Estimates 


Term  Estimate  Sid  Error  t  Ratio  Prob>jtj 

Intercept  0.0262436  0.00105  24.09  <,0001 

Proto{R2|  0.006568?  0.001461  4.50  <.0001 

ProiotRSl  -0.060111  0.001484  -0.0?  0.0410  _ 

[  Effect  Tests  .  . ^ . . 

Source  Nparm  OF  Sum  of  Squares  F  Ratio  Pmb>F 

Proto  2  2  0.00127087  12.8360  <.0001 
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Appendix  F  (Continued) 


Rej  Response  to  Protocol  (Metabolic  Demand)  by  Ensemble 


fEnsemble^A 


Response  ReT 


Observations  (or Sam  Vgjts)  44 

[  Analysis  of  Variance 


Source  Df  Sum  ot  Squares  Mean  Square  F  Ratio 


Model  a  0.00002880  0000014  28115 

Error  41  000021000  0.000005  Prob>F 

C.  Total  43  0.00023880  _ _  0.0717 

[Parameter  Estimates  I 


Term  Estimate  Sid  Error  l  Ratio  Rrob>lf| 

Inter^pl  0.0114238  0.000341  33.47  <0001 

M|M1|  *0.000352  0, 0004 8B  0.72  0.4740 

#§  _  0.001109$  0.00048  2.31  0.0259 

(Effect Tests  “  '  '  T  _  ) 

Source  Npami  DF  Sam  of  Squares  F  Ratio  Ffot>>F 

M  2  2  0.00002880  28115  0.0717 


Appendix  F  (Continued) 


Re,T  Response  to  Protocol  (Metabolic  Demand)  by  Ensemble 


Root  Mean  Square.  Error  0,00271  5 
Mean  of  Response  0,01 2167 

Observations  (or  Sum  V^ts)  42 


M3 


0,01114286  0,00072584  0.011143 


[Analysis  of  Variance 


Source  DF  Sum  of  Squares  Moan  Square  F  Ratio 

Model  2  O.OOO04433  O.OOO022  3.8070 

Error  39  0.80028758  0.000007  Pret>>F 

C,  Total  41  0.00033183  0  8810 


Parameter  Estimates 


Term  Estimate  Std  Error  I  Ratio  Pfob>|tj 

Intercept  0.012168?  0.000419  29.04  <8001 

Mpll  0.0014048  0.000892  2.37  0.0228 

IM  *0  >000381  0.00D592  43.64  0.5240 

[Effect  Tests 


Source 

M 


Nparm 

2 


OF  Sum  of  Squares 
2  0.0O0O4433 


F  Ratio  Prob  >  F 
3.0070  0.0610 
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Appendix  F  (Continued) 


Re/r  Response  to  Protocol  (Metabolic  Demand)  by  Ensemble 


Observations  (or  Sum  VVpls)  40 


( Analysis  of  Variance  j 

Source 

DP  Sum  of  Squares 

Mean  Square  P  Patio 

Model 

2  0.00013011 

0.000070  10.7741 

Error 

43  0,00027780 

0.000006  Prob  >  p 

C.  Total 

45  0.00041872 

O.0OO2 

{ Parameter  Estimates 

j 

Term 

Estimate  Std  Error 

matte 

Pmb>Ji| 

Intercept 

0,0125782  0.000375 

33,56 

<,0001 

MfMll 

0.0023583  0,000524 

4.30 

<,0001 

MIM2) 

-0.000846  0.000533 

*1 ,21 

0.2321 

f  Effect  Tests  j 

Source  Nparm  DF  Sum  of  Squares  f  Ratio  Prob  >  F 
M  2t  0.00013811  10.7741  0.0002 


140 


Appendix  F  (Continued) 


Re,T  Response  to  Protocol  (Metabolic  Demand)  by  Ensemble 


( £n$emb!e~D 

j  Response  ReT 

{ Whole  Model 

Tii/i' . 

A_  . . . . . . . 

- — - — . k.skksssc^ 

[  Actual  by  Predicted  Plot 

[Leverage  Plot  1 

-i — i — i — i — . — | — . — | — f 

.010  .015  .020  .025  .030  .035 1 

Ret  Predicted  P<.0001  R$<p0.36 
RMSE*0.0O33 


0.035- 


g  0.030- 
!  0,025«| 


0.020-4 


'O.OfO-K 


f  Summary  of  Fit 

RSquare 
RSquare  Adj 
Root  Mean  Square  Error 
Mean  of  Response 


012  .013  .014  .015  ,016  .017  .018  .0191 
M  leverage,  P<  0001 

j  toast  Squares  Moans  Table 


0.364402 

0,33423 

0.003349 

0,015280 

45 


Level  least  Sq  Mean 

m  0.01833333 

M2  0.01520000 

M3  0.01233333 


Std  Error  Mean 
0.00080471  0.018333 
0.00080471  0.015200 
0.00088471  0,012333 


[Analysis 

of  Variance 

: ) 

Source 

DF  Sum  of  Squares 

Mean  Square  f  Ratio 

Model 

2  0.00027018 

0.000135  12.0444 

Error 

42  0.0004710? 

0.000011  Preb  >  F 

C.  Total 

44  O.OO074124 

<.0001 

|  Parameter  Estimates 

Term 

Estimate  Std  Error 

1  Ratio 

Prob>|t| 

Intercept 

0.0152889  0.000490 

30.62 

<.0001 

MfMl| 

0.0030 444  0.000706 

4.31 

<0001 

MJM2J 

-0.000089  0,000706 

-0.13 

0.9004 

(  Effect  Tests  1 

Source  Nparm  DF  Sum  of  Squares  F  Ratio  Prob>F 

M  2  2  0.00027018  12.0444  <.0001 
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Appendix  F  (Continued) 


Re/r  Response  to  Protocol  (Metabolic  Demand)  by  Ensemble 


Observations  (or  Sum Wgts)  48 


{ Analysis  of  Variance 

1 

Source  DF  Sum  of  Squares  Mean  Square  ! 

Model  2  0.00072507  0.000383  2 

Error  4$  0.0007SS70  0,000017  R 

C,  Total  47  0.00148187 

f  Ratio 
18148 
rob  >  F 
<.0001 

! 

I  Parameter  Estimates 

Term  Estimate  SW  Error  t  Ratio  Pmb>;i| 

Intercept  0.0236888  0.000592  40.00  <0001 

M(M1)  0,0045614  0.000837  5.45  <.0001 

MJM21  0.0002448  0.000851  0.29  0.7749 

|  Effect  Tests  j 

Source  Nparrrt  DF  Sum  of  Squares  F  Ratio  Rrob>F 

M  2  2  O.0QO725&7  218146  <.0001 
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APPENDIX  G 

JMP  IN  DATA  ANALYSIS  -  M2R5 
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Appendix  G 

Rej  Response  to  Ensemble  by  Phase  (M2R5  Dataset) 


Mean  of  Response  0.0128$ 

Observations  (or  Sum  Wgts)  20 

f  Analysis  of  Variance  ; ‘ ;  , 

Source  OF  Sum  of  Squares  Mean  Square  FftatJo 

Model  1  000000076  7.502807  0.0878 

Error  27  0.00023S45  0,00O68$8  Frob>F 

C.  Total  28  0.00023421  _  0.7602 

Parameter  Estimates _ ] 

Term  Estimate  Std  Error  t  Ratio  FroO>|t| 

Incept  0.0126952  0,000548  28.24  <M0i 

YrfTll  O.OO0<815  _  0.000546  0.80  0.7602 

[effect  Tests  . . .  ~ .  3 

Source  Npsrm  OF  Sum  or  Scares  F  Ratio  Frot>  >  F 

Vr  1  1  7.592780-7  00878  0.7692 
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Appendix  G  (Continued) 


R^t  Response  to  Ensemble  by  Phase  (M2R5  Dataset) 


RSquare  Adj  -&0t718  Y2  0.01178571  0.0008B665  0.011785 

Root  Mean  Square  Erro r  0.00301 0 

Mean  of  Response  0.01220? 

Observations  (or  Sum  Vtffis) _ 29  _ 

[  Analysis  of  Variance  ~  j 


Source  Of  Sum  of  Squares  Mean  Square  f  Ratio 

Model  1  0.00000430  0.000004$  0.5271 

Emt  2?  0.000245 06  0.6000091  Pm  >  F 

C,  Total . jj»  0.00025076  04741 

[iparamater  Estimates  "" . : ""  •] 

Term  Estimate  Std  Error  t  Ratio  Frob>|t| 

Intercept  0.0121029  0.000561  21.74  <0001 

YifyiJ  0.0004071  0.000561  0.73  0.4741 


Effect  Tests 


□ 


Source  Nparm  DF  Sum  of  Squares  F  Ratio  Pret)>F 

Vf  1  1  0.0000048  0.5271  0,4741 
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Appendix  G  (Continued) 


Re}T  Response  to  Ensemble  by  Phase  (M2R5  Dataset) 


( Ensembie^C 


I  Response  Re? 


1  Whole  Model 


Yr 


[Aiuallbyiredlcted  Riot 


[  Leverage  Riot 


0.030- 


Re?  Predicted  P*mm  RSqs0.15 

RMSE^0.0035 

[Summary  of  Fit  '  '  j 

R  Square  0.154S47 

PSquamAdj  0.123234 

Root  mm  Squam  Error  0.003542 

mm  of  Response  0.01334$ 

Observations  (or  Sum  Wgts)  2$  _____ 

[  Analysis  of  Variance  ' '  . . "  . ; 

Source  OF  Sum  of  Squares  Mean  Square  F  Ratio 

Mode!  1  0 .000081 83  0.000062  4.035$ 

Error  27  0,000330$$  0.000013  Rrob>E 

C,  Total  2$  0.000400$$  0.0340 


J  (MB* 


#  r  * 


— i — i — i — i — i — i — I 

.0115,012.01 2$  ,013.0135  .014,0145.015 

Yr  Leverage.  Ps0.0349 

[  toast  Squares  Means  Table  _ •; . JT~ 

Level  LcastSq  Mean  Sid  Eiror  Mean 
Y1  0.01485714  0.000046 52  0.014857 

Y2  0.01103333  0.00001442  0.011933 


[parameter  Estimates 


Teem  Estimate  Std  Error  I  Ratio  Prob>ll| 

Intercept  00133952  0.0906$$  20.36  <.0001 

YrfYlf  0.0014619  0.0OC658  2.22  0.034$ 


f  Effect  Tests 


Source 

Yr 


Nparm 

1 


OF  Sum  of  Seisms 
1  0,00006100 


F  Ratio  Prob>F 
4.9355  0.0649 
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Appendix  G  (Continued) 


APPENDIX  H 

GRAPHS  OF  Re, t  VERSUS  AP  BY  ENSEMBLE 
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Appendix  H  (Continued) 


Graphs  of  R^t  versus  AP  by  Ensemble 


Appendix  H  (Continued) 


Graphs  of  R^j  versus  AP  by  Ensemble 


Appendix  H  (Continued) 
Graphs  of  Re, t  versus  AP  by  Ensemble 

Phase  1  Ensemble  E 
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